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ABSTRACT 

Context. The massive and luminous star-forming region W49A is a well known Galactic candidate to probe the physical conditions 
and chemistry similar to those expected in external starburst galaxies. 

Aims. We aim to probe the physical and chemical structure of W49A on a spatial scale of ~0.8 pc based on the JCMT Spectral Legacy 
Survey, which covers the frequency range between 330 and 373 GHz. 

Methods. The wide 2x2 arcminutes field and the high spectral resolution of the HARP instrument on JCMT provides information 
on the spatial structure and kinematics of the cloud traced by the observed molecular lines. For species where multiple transitions are 
available, we estimate excitation temperatures and column densities using a population diagram method that takes beam dilution and 
optical depth corrections into account. 

Results. We detected 255 transitions corresponding to 60 species in the 330-373 GHz range at the center position of W49A. Excitation 
conditions can be probed for 16 molecules, including the complex organic molecules CH 3 CCH, CH 3 CN, and CH 3 OH. The chemical 
composition suggests the importance of shock-, PDR-, and hot core chemistry. Many molecular lines show a significant spatial extent 
across the maps including CO and its isotopologues, high density tracers (e.g. HCN, HNC, CS, HCO^), and tracers of UV-irradiation 
(e.g. CN and C 2 H). The spatially extended species reveal a complex velocity-structure of W49A with possible infall and outflow 
motions. Large variations are seen between the sub-regions with mostly blue-shifted emission toward the Eastern tail, mostly red- 
shifted emission toward the Northern clump, and emission peaking around the expected source velocity toward the South-west clump. 
Conclusions. A comparison of column density ratios of characteristic species observed toward W49A to Galactic PDRs suggests 
that while the chemistry toward the W49A center is driven by a combination of UV-irradiation and shocks, UV-irradiation dominates 
for the Northern Clump, Eastern tail, and South-west clump regions. A comparison to a starburst galaxy and an AGN suggests similar 
C 2 H, CN, and H 2 CO abundances (with respect to the dense gas tracer ^^'CS) between the ~0.8 pc scale probed for W49A and the >1 
kpc regions in external galaxies with global star-formation. 
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1. Introduction 

Star-formation in galaxies occurs on different scales from iso¬ 
lated cores / globules with sizes of <1 pc through giant molec¬ 
ular clouds of a few 10 pc up to the scales of starburst galaxies 
representing star-formation on scales of a few 100 pc and with 
luminosities of Gr > lO" L© for luminous infrared galaxies 
(LIRGs) and Ljr > 10^^ L© for ultraluminous infrared galax¬ 
ies (ULIRGs). These regions also represent a large spread in 
star-formation rates with the 3 Mo/yr in the Milky Way to ~10^ 
Mo/yr in starburst galaxies (e.g. Solomon & Vanden Bout 2005). 
The different types of star-forming regions in galaxies may trace 
different physical processes that control the star-formation in the 
different environments. 

The Galactic starburst-analogue W49A with its more than 
~ 100 pc size (Simon et al. 2001) represents an intermediate case 
between the star-formation seen in Galactic star-forming regions 
and in starburst galaxies. With its intense star-formation com¬ 


pared to other Galactic sources, studying its physical and chem¬ 
ical properties can give insights into the mechanism behind the 
starburst phenomenon. In our previous work based on part of the 
data presented in this paper we focused on the physical proper¬ 
ties of W49A (Nagy et al. 2012). In this work we focus more on 
the chemistry of this Galactic starburst analogue. W49A is one 
of the most massive (M ~ 10® M©, Sievers et al. 1991) and lumi¬ 
nous (>10^ L©, Ward-Thompson & Robson 1990) star-forming 
regions in the Galaxy. Even though its luminosity is lower than 
that of starburst galaxies, W49A is more comparable to ULIRGs 
in terms of its luminosity per unit mass of gas. In ULIRGs values 
of ~100 Lq / M© are measured (e.g. Solomon et al. 1997) while 
the value for W49A is >10 Lq / M©. The luminosity of W49A 
is a result of an embedded stellar cluster containing the equiva¬ 
lent of about 100 07 stars (Conti & Blum 2002), corresponding 
to four stellar clusters in a 5'x5' (16x16 pc) region around the 
center based on J, H, and Ks images (Alves & Homeier 2003). 
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Wu et al. (2014) discovered a very massive star corresponding 
to the central cluster of W49A, and estimated its mass to be in 
the range between 90 and 250 Mq. Apart from the stellar pop¬ 
ulation observed at near-infrared wavelengths, the ongoing star 
formation which appears to be forming a comparable mass of 
stars can be studied in the radio continuum (e.g. De Free et al. 
1997, 2004, 2005) and at mm- and sub-mm (e.g. Wilner et al. 
2001, Galvan-Madrid et al. 2013) wavelengths as hot cores and 
Ultra-Compact (UC) Hit regions. 

The distance of W49A was recently estimated to be 
ll.lUg™ kpc (Zhang et al. 2013), which is consistent with the 
11.4 kpc value derived by Gwinn et al. (1992). Due to the large 
distance of W49A its small-scale structure including hot cores, 
outflows, and UC Hii regions require interferometric observa¬ 
tions to be fully resolved (e.g. Wilner et al. 2001, De Free et al. 
2005, Galvan-Madrid et al. 2013). However, previous studies us¬ 
ing single-dish telescopes have already provided information on 
the chemical complexity of W49A, and the physical and chemi¬ 
cal properties of the region. 

HCO^ 1-0 lines with red-shifted absorption and blue-shifted 
emission were observed by Welch et al. (1987) and interpreted 
as evidence of global collapse toward the central ~2 pc region 
of W49A. An alternative explanation for the large number of 
O-stars contributing to the luminosity of W49A was proposed 
by Serabyn et al. (1993), who suggested a cloud-cloud collision 
based on multiple transitions of CS and C^^S (from 7=3-2 to 
7=10-9). However, based on more recent observations, this ex¬ 
planation is less probable than the global cloud collapse sce¬ 
nario. Large-scale maps of W49A (e.g. Galvan-Madrid et al. 
2013) show a hierarchical network of filaments, converging from 
larger scales to the center of W49A (W49N). 

Signatures of the high star-formation activity were discov¬ 
ered using Spitzer mid-IR images and position-velocity dia¬ 
grams based on *^CO 7 = 2-1 and C'^O 7 = 2-1 maps 
by Feng et al. (2010), who have identified two expanding shells 
toward W49A. 

Vastel et al. (2001) have studied the physical conditions in 
the FDR component of W49A using the FIR lines of [Oi] and 
[Cii] observed with the Long Wavelength Spectrometer of the 
Infrared Space Observatory as well as rotational lines of CO (CO 
1-0, C'*0 2-1) observed with the 15-m SEST telescope, and de¬ 
rived a radiation held of Go = 3 x 10^ (in units of 1.3x10“^ erg 
s“' cm“^ sr“') and an average gas density of 10 "^ cm“^. 

A recent study based on data from the JCMT Spectral 
Legacy Survey (SLS, Flume et al. 2007) focused on an extended 
warm and dense gas component toward W49A seen in H 2 CO 
(Nagy et al. 2012), and characterized the physical properties and 
excitation of the region using CH 3 OH, SO 2 , H 2 CO, and HCN 
transitions. Another study that makes use of the SLS line sur¬ 
vey data focuses on dense gas tracers (Roberts et al. 2011), such 
as HCN, HNC, DCN, HCO^, and their isotopologues, and by 
a comparison to line ratios measured toward starburst galaxies 
and Active Galactic Nuclei (AGN), as well as to chemical mod¬ 
els, finds that W49A is a template for starburst galaxies rather 
than for AGN. 

In this paper, we present results on the chemical inven¬ 
tory of W49A based on the SLS, carried out with the James 
Clerk Maxwell Telescope (JCMT) at a resolution of ~15". This 
line survey provides the largest frequency coverage data toward 
W49A up to date, and as such, it provides a useful starting point 
for future higher resolution studies with instruments such as 
ALMA. 


2. Observations and data reduction 

The SLS observations have been carried out using the 16- 
receptor (spatial pixel) Heterodyne Array Receiver Frogramme 
B (HARF-B, 325-375 GHz) and the Auto-Correlation Spectral 
Imaging System (ACSIS) correlator (Buckle et al. 2009) at 
the James Clerk Maxwell Telescope^ (JCMT) on Mauna Kea, 
Hawai’i. The observations were carried out in jiggle position 
switch mode, sampled every 7.5" for a 2x2 arcminutes field cen¬ 
tered on RA(J2000) = 19‘'10™13.M; Dec(J2000) = 09°06'14". 
The spectra were calibrated using an off-position correspond¬ 
ing to an offset (-h 840", -1-840") compared to the central posi¬ 
tion. Observations toward the reference position have been car¬ 
ried out to confirm that the off position is reliable and is without 
emission at the observed frequencies. The pointing was checked 
every hour and is estimated to be accurate to 1.5". The angu¬ 
lar resolution of the JCMT is ~15" at 345 GHz, equivalent to 
~0.8 pc at the distance of W49. The spectral resolution at the 
observed frequencies is ~ 0.8 km s ' and the beam efficiency is 
0.63 (Buckle et al. 2009). The original line survey was carried 
out in the 330-360 GHz frequency range and was later extended 
to 373 GHz. This paper summarizes the results from the whole 
330-373 GHz line survey. 

The data reduction was done using a combination of tasks 
from the Starlink package and the ORAC Data Reduction 
pipeline (ORAC-DR). We reduced the observed time-series 
cubes using the ORAC-DR pipeline, which creates three- 
dimensional cubes after checking for consistency between the 
calculated Tsys and observed rms noise in the data; checking for 
variations in the rms noise measured by each receptor across the 
map, removes baselines from every spectrum and co-adds spec¬ 
tra that corresponds to the same position and frequency. The re¬ 
sults of the pipeline were checked and were corrected for re¬ 
maining bad data, such as spectra with high rms noise level, bad 
baselines, and spectral ranges affected by spikes, using a combi¬ 
nation of Starlink tasks. 

Figure 1 shows the typical rms noise levels as a function of 
frequency in the 330-373 GHz frequency range measured at ev¬ 
ery 2 GHz. Typical noise levels in units are in the range 
between 0.02 and 0.12 K. The highest noise level measured at 
~368 GHz is due to the poor atmospheric transmission as shown 
in the top panel of Fig. 1. Atmospheric transmissions are shown 
for typical precipitable water vapor (FWV) levels of 1 mm and 
3.75 mm on Mauna Kea. 

3. Results 

3.1. Line identification 

We have identified the lines detected toward the center of W49A 
(RA(J2000) = 19’'10”13.M; Dec(J2000) = 09°06'14") in the 
frequency range between 330 and 373 GHz. Analysed posi¬ 
tions other then the center include the eastern tail (RA(J2000) 
= lOMO^lb.^b; Dec(J2000) = 09°05'48"), northern clump 
(RA(J2000) = 1 o'* 10” 13."6; Dec(J2000) = 09°06'48") and 
southwest clump (RA(J2000) = 19''10”10."6; Dec(J2000) = 
09°05'18"). These positions were selected based on differ¬ 
ent kinematical signatures traced by dense gas tracers such 
as HCO+, HCN, and HNC (Roberts et al. 2011) and are 
shown for example on Fig. 3. Appendix A shows details on 

' The James Clerk Maxwell Telescope is operated by the Joint 
Astronomy Centre on behalf of the Science and Technology Facilities 
Council of the United Kingdom, the Netherlands Organisation for 
Scientific Research, and the National Research Council of Canada. 
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Fig. 1. The typical rms noise levels in T* units measured in the 
SLS survey as a function of the frequency (bottom panel) and 
the corresponding atmospheric transmission for typical levels of 
water vapor at Mauna Kea (top panel). 


resolved blends and have contribution to the line emission from 
other molecules / transitions. However the excitation diagram 
analysis provides a cross-check on the line identifications (and 
indeed highlighted some initial mis-assignments) and we believe 
that misindentified lines are not a signifiant issue. 

To measure the spatial extent of the detected species, we ap¬ 
ply a 2-dimensional Gaussian fit to the integrated line intensi¬ 
ties, using the BEAMFIT task in Starlink. The values shown in 
Table 1 correspond to the FWHM values of the fit. For molecules 
with multiple detected transitions we include their largest spatial 
extent. For molecules with measured source sizes of ~23x23", 
the lines were only detected toward a few positions around 
the source center. In the following sections we quote sizes of 
~ 17x17" for these species, to correct for the ~15" beam size 
of JCMT around the observed frequencies. The center positions 
of the Gaussian fits corresponding to the different species are 
within a few arcsecond distance around the W49A center. In 
addition to CO and its isotopologues, high density tracers (e.g. 
HCN, HNC, CS, HCO+) and tracers of UV-irradiation (e.g. CN 
and C 2 H) show the largest spatial extent (Table 1). The mea¬ 
sured position angles cover a large range indicating a complex 
source structure. Section 3.2 includes a summary of the detected 
molecules. We discuss the excitation conditions and column den¬ 
sities of the species with multiple detected transitions in Section 
3.4. 


the detected species. We added a comment when the lines 
were detected in the high-mass protostar AFGF 2591, another 
source in the JCMT SFS, based on Van der Wiel (2011). 
The spectroscopic data are based on the Cologne Database for 
Molecular Spectroscopy (CDMS, Muller et al. 2005)^ and the 
Jet Propulsion Faboratory (JPF, Pickett et al. 1998)^ molecular 
databases. The identification is based on an initial search range 
of transitions with an upper level energy of 400 K, and was ex¬ 
tended in particular cases such as vibrationally excited lines of 
HCN and HNC. In the case of weak (<3cr detections) lines, an 
additional check of the spatial distribution and a cross-check of 
the detection in data corresponding to different rest frequencies 
was applied to avoid possible artifacts e.g. spurs or baseline er¬ 
rors. For the identification we considered velocities in a range 
between +5 km s“'around the expected source velocity of 10 
km s 'and extended it for the most asymmetric lines. 

Toward the centre position we detected 255 lines in the 330- 
373 GHz frequency range, that belong to 60 molecular species 
summarized in Table 1. Most of the detected lines are related 
to sulphur-bearing species, such as SO 2 and its isotopologues. 
Various species are detected that indicate the possible impor¬ 
tance of shock- and hot core chemistry and UV irradiation (see 
Sect. 4). The detected transitions cover a large energy range up 
to ~1067 K corresponding to vibrationally excited HCN (Fuller 
et al, in prep.). However, most transitions have upper level ener¬ 
gies below 400 K (Fig. 2). Excitation conditions can be probed 
for 16 species (Table 2) based on multiple detected transitions 
with a sufficient energy range. Among the 255 detected lines, 
there are 3 unidentified lines (U-lines in Appendix A.l). With 
the large number of spectral features misassignment to species 
is possible. Also some of the line profiles are non-Gaussian, and 
the FWHM (full width at half maximum) width of the lines vary 
over a relatively large range (~7-19 km s“'toward the central po¬ 
sition). This may suggest that some of the identified lines are un- 

^ http: //WWW. astro. uni-koeln. de/cdms/catalog 

^ http://spec.jpl.nasa.gov 



Fig. 2. The number of transitions detected versus the upper state 
energy of the transition in 25 K wide bins. 


3.2. The detected species 

The detected species summarized in Table 1 can be related to 
various processes including shocks, UV-irradiation by the em¬ 
bedded OB stars, and to hot core chemistry that includes com¬ 
plex organic molecules that are released from the grains into the 
gas-phase at temperatures of 100-300 K and are shielded from 
the dissociative UV-radiation. 

3.2.1. Shock tracers 

A large number of the detected molecules can be related to shock 
chemistry (Fig. 3). All of these molecules are observed toward a 
< 20" X 20" field around the center, but not detected toward the 
other subregions analysed in this paper. 
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Table 1. The measured spatial extent of the species detected toward the center of W49A. The position angle is measured from North 
through East. 


Species 

Number of 
transitions 

£i,p 

(K) 

Major 
axis (") 

Spatial extent® 

Minor Pos. angle'’ 

axis(") (deg) 

Species 

Number of 
transitions 

£up 

(K) 

Major 
axis(") 

Spatial extent® 
Minor Pos 

axis (") 

angle 

(deg) 

CO 

1 

33 

98 

60 

90 

"CO 

1 

32 

72 

54 

87 

c‘’o 

1 

32 

39 

28 

69 

CjH 

3 

42 

61 

46 

79 

CN 

3 

33 

54 

39 

77 

HCN 

1 

43 

37 

27 

76 

cs 

1 

66 

39 

24 

70 

DCN 

1 

52 

36 

23 

146 

CH 30 H 

26 

17 - 372 

30 

23 

64 

H 2 CO 

7 

52 - 241 

29 

23 

64 


1 

65 

28 

23 

69 

H 2 CS 

6 

91 - 209 

30 

23 

142 

NS 

1 

70 

29 

23 

9 

NO 

3 

36 - 209 

29 

23 

62 

^^SO 

3 

78-87 

33 

23 

112 

HNC 

1 

43 

31 

23 

67 

HCO+ 

1 

43 

32 

26 

179 

H 

2 


35 

23 

77 

SO 2 

41 

31-612 

23 

23 


N 2 H+ 

1 

45 

23 

23 


^-’S 02 

18 

35 - 338 

23 

23 


‘"CS 

1 

80 

23 

23 


^''S 02 

33 

35 - 547 

23 

23 


H30+ 

1 

140 

23 

23 


S0‘®0 

4 

117-326 

23 

23 


SiO 

1 

75 

23 

23 


so"o 

1 

58-180 

23 

23 


SO 2 , V 2 =l 

9 

805 - 998 

23 

23 



1 

45 

23 

23 


SO 

6 

26 - 143 

23 

23 


HCS* 

1 

74 

23 

23 


"“so 

5 

25-85 

23 

23 


H'^CN 

1 

41 

23 

23 


OCS 

3 

237-271 

25 

23 


SO+ 

1 

70 

23 

23 


HNCO 

7 

127 - 204 

23 

23 


0 ‘^C^'‘S 

1 

247 

23 

23 


CH 3 CN 

5 

151 -215 

23 

23 


HC'^N 

1 

41 

23 

23 


H‘"CO"' 

1 

42 

23 

23 


HN'^C 

1 

42 

23 

23 


CH 3 CCH 

4 

172 - 254 

23 

23 


HCN, V 2 =l 

1 

1067 

23 

23 


HC 3 N 

3 

307 - 376 

23 

23 


HC 3 N, V 7 = 1 

4 

629-663 

23 

23 


HCO 

1 

74 

23 

23 


CH 3 OCH 3 

2 

48-49 

23 

23 


s‘®o 

3 

91-99 

23 

23 


CO+ 

1 

34 

23 

23 


H 2 S 

2 

135 - 263 

23 

23 


S"0 

2 

76 - 103 

23 

23 


H("CO 

2 

61-65 

23 

23 


CH 3 CHO 

1 

155 

23 

23 


HNC, V 2 =l 

2 

709-710 

23 

23 


H'^NC 

1 

43 

23 

23 


HC‘’ 0 + 

1 

42 

23 

23 


H 2 CN 

3 

53 - 100 

23 

23 


HO‘"C+ 

1 

41 

23 

23 


HC‘*0+ 

1 

41 

23 

23 


H 2 C"“S 

1 

118 

23 

23 


CH 3 COCH 3 

1 

282 

23 

23 


Hj"S 

1 

154 

23 

23 


HfS 

1 

154 

23 

23 









“ The error corresponding to the fitted major and minor axis values is a few arcseconds. 

’’ Position angles are only derived for species with significant spatial extent and elongation in order to minimise the uncertainty. 


- SiO was detected in hot and shocked regions, such as molec¬ 
ular outflows (e.g. Nisini et al. 2007) and supernova rem¬ 
nants (e.g. Ziurys et al. 1989). The observations suggest that 
silicon is released from grain mantles in these regions. A pro¬ 
cess that can release silicon into the gas phase is the sput¬ 
tering of (charged) grains by heavy neutral particles in C- 
shocks (Schilke et al. 1997, Gusdorf et al. 2008). SiO is de¬ 
tected toward W49A with a spatial extent of ~ 17" x 17" (cor¬ 
rected for the 15" beam size of JCMT) covering the positions 
corresponding to the central stellar cluster in its 7 = 8 - 7 
transition. Its 7 = 2 - 1 transition has previously been de¬ 
tected toward W49A by Lucas & Liszt (2000). 

- Most of the detected species are sulphur-bearing molecules, 
such as SO and SO 2 and their isotopologues ^^S02, ^^S02, 
SO'^O, SO'^O, ^‘^SO, 33sO, S^^O, and S'^O. Sulphur is 
frozen onto grain mantles and can be sputtered in shocks, 
leading to the formation of the species mentioned above (e.g. 
Wakelam et al. 2004). We also detected 8 transitions of vi- 
brationally excited SO 2 . These molecules are only detected 
toward the center of W49A and spatially confined to a region 
of - 17" X 17", similar to SiO. 


- Other detected sulphur-bearing species include H 2 S, OCS 
(and its isotopologue 0'^C^"^S). H 2 S and OCS show spatial 
extents similar to that of SO 2 , SO, and their isotopologues; 
and to SiO. 

- We also detected SO^, which has been reported to be a diag¬ 
nostic of dissociative shock chemistry (Turner 1992), but has 
also been detected toward various PDRs (see Sect. 3.2.2). 

- We detected 6 transitions of HNCO, which has also been 
reported to trace shocks toward the Galactic Center and in 
molecular outflows, such as LI 157 (Rodriguez-Lernandez et 
al. 2010). 
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Fig. 3. The spatial distribution of species that can be related to shock chemistry. The crosses correspond to the subregions discussed 
in this paper: source center, eastern tail, northern clump, and south-west clump. The contours correspond to the 20%, 40%, 60%, 
and 80% of the maximum line intensity. The SO 2 , V 2 =l integrated intensity was calculated as the total integrated intensity of the 
three blended SO 2 , V 2 =l lines around 366.1 GHz. 
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Fig. 4. The spatial distribution of species that can be related to PDR chemistry. The crosses correspond to the subregions discussed 
in this paper: source center, eastern tail, northern clump, and south-west clump. For C 2 H and CN the contours correspond to the 
30%, 40%, 60%, and 80% of the maximum line intensity. For CO^ and HCO the contours correspond to the 45%, 60%, and 80% 
of the maximum line intensity. The different percentages used for the contour levels for the CO^ and HCO integrated intensities is 
due to the lower signal-to-noise of these species compared to C 2 H and CN. 
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3.2.2. Tracers of UV-irradiation 

The embedded clusters of O and B stars in the W49 region (e.g. 
Alves & Homeier 2003) are expected to create PDRs. Several 
species that can be related to PDR chemistry are detected in the 
SLS frequency range (Fig. 4), some of them with a significant 
spatial extent. 

- We detected CN with a large extent (~54" x 39"), includ¬ 
ing all subregions. This suggests the possible importance of 
PDR chemistry in W49A and may be tested using the ratio of 
CN/HCN, which is a well-known tracer of PDRs (e.g. Fuente 
et al. 1996). The spatial distribution of HCN and its iso- 
topomer HNC based on SLS data were analysed in Roberts 
et al. ( 2011 ). 

- C 2 H (ethynyl) is a commonly observed molecule in the inter¬ 
stellar medium, including in PDRs. In star-forming regions 
it was proposed to be related to the earliest stages of mas¬ 
sive star-formation (Beuther et al., 2008). We detected three 
transitions toward W49A, with a very large spatial extent 
(~61" x46"). 

- HCO (formyl radical) is also detected in a ~17" x 17" re¬ 
gion around the center. HCO was found to be a tracer of illu¬ 
minated cloud interfaces as it was detected in several PDRs 
including the Orion Bar, NGC 2023, NGC 7023, and S140 
(Schilke et al. 2001) and the Horsehead PDR (Gerin et al. 
2009). 

- The SO^ and CO^ reactive ions have been detected in var¬ 
ious PDRs, such as the Orion Bar (Fuente et al. 2003), and 
in other regions of high FUV- and X-ray irradiation. We de¬ 
tected both CO^ and SO^ toward the center of W49A with a 
spatial extent of ~17" x 17". 

3.2.3. Complex organic molecules 

We have detected a number of complex organic molecules (e.g. 
Herbst & Van Dishoeck 2009), that are most likely related to hot 
cores. 

- CH 3 OH is detected over a 30" x 20" region around the cen¬ 
ter. Its spatial distribution is shown in Nagy et al. (2012). Its 
'^CH 30 H isotopologue is also detected. 

- We detected 5 transitions of CH 3 CN (methyl cyanide), that 
has been detected in hot cores, including Orion KL. The 
7 = 12 - 11 transition of CH 3 CN was previously detected 
toward W49 using the BIMA array and is likely connected 
to hot cores and also follows the compact dust concentrations 
(Wilneretal. 2001). 

- We detected four transitions of HC 3 N (cyanoacetylene), that 
has been detected toward various Galactic star forming re¬ 
gions as well as in a number of dust-rich galaxies (Lindberg 
et al. 2011). As it is destroyed by UV radiation and in re¬ 
actions with the ion, HC 3 N is expected to trace warm, 
dense, and shielded regions. This may explain the fact that 
HC 3 N is seen to be confined toward the central ~17" x 17" 
region. Some of the vibrational transitions of HC 3 N (vy^l) 
were also detected toward the W49A center. Vibrationally 
excited HC 3 N has also been detected toward several hot 
cores (Wyrowski et al. 1999), as well as the galaxy NGC4418 
(Costagliola & Aalto 2010). 

- One transition of CH 3 CHO (acetaldehyde) are detected to¬ 
ward the W49A center. A previous survey by Ikeda et al. 
( 2001 ) detected this molecule toward various hot molecular 
clouds. 


- Four transitions of CH 3 CCH (propyne) are detected in our 
line survey at a few central pixels with a spatial extent of 
~17" X 17". This molecule has been detected in a number of 
starburst galaxies, such as M82 (Aladro et al. 2011a). 



40 20 0 -20 -40 40 20 0 -20 -40 

ARA ["] 


Fig. 5. The spatial distribution of the detected most spatially ex¬ 
tended complex organic molecules. The crosses correspond to 
the subregions discussed in this paper: source center, eastern tail, 
northern clump, and south-west clump. The contours correspond 
to the 20%, 40%, 60%, and 80% of the maximum line intensity. 


3.2.4. Other detected species 

- CO, '^CO, and C'^O are covered by the SLS survey, and are 
detected at each position of the SLS field (Figure 7), allowing 
us to probe the gas kinematics and column density. 

- The only deuterated species we detected toward W49A is 
DCN (Roberts et al. 2011). An upper limit for the covered 
but not detected DCO^ 7=5-4 transition was also reported by 
Roberts et al. (2011). A low abundance of deuterated ions is 
expected as deuterated molecular ions are destroyed rapidly 
at gas temperatures above 30 K (e.g. Roberts & Millar 2007). 

- H 2 CO (formaldehyde) is a commonly used temperature and 
density tracer in Galactic star-forming regions (e.g. Mangum 
& Wootten 1993) and in external galaxies (e.g. Mangum et 
al. 2008) as well. Previous results from this line survey fo¬ 
cused on this molecule, which revealed a large ~3x3 pc re¬ 
gion around the W49A center with kinetic temperatures of 
>100 K and densities of >10^ cm“^ (Nagy et al. 2012). Its 
isotopologue Hj^CO is also detected toward the center. 

- Similar to H 2 CO, H 2 CS (thio-formaldehyde) is a near- 
prolate rotor, which allows estimation of both kinetic tem¬ 
perature and density. Though its spatial extent (30" x 17") 
given by a 2D Gaussian fit toward the center is similar to that 
of H 2 CO, H 2 CS was not detected toward the other main re¬ 
gions (Eastern tail. Northern clump, and Southwest clump). 
Its isotopologue H 2 C^"^S is also detected. 

- We detected a few ions, including HCO^ and its isotopo- 
logues HC'^O^, H^^CO^, and HC^^O^. As a dense gas 
tracer, HCO^ data in the SLS survey were analysed by 
Roberts et al. (2011). Other detected ions not mentioned in 
the sections above include H 30 ^ and N 2 H^. N 2 H^ shows 
a larger spatial extent than H 30 ^ covering the Eastern tail 
and Northern clump regions. H 30 ^ is an important chemi¬ 
cal ingredient of molecular clouds as it can be related to the 
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Fig. 6. The spatial distribution of N 2 H^ and H 30 ^ ions and CS. The crosses correspond to the subregions discussed in this paper: 
source center, eastern tail (ET), northern clump (NC), and south-west clump (SW). The contours correspond to the 20%, 40%, 60%, 
and 80% of the maximum line intensity. 



ARA ["] 


Fig. 7. The spatial distribution of CO, *^CO, and line emission. The crosses correspond to the subregions discussed in this 
paper: source center, eastern tail (ET), northern clump (NC), and south-west clump (SW). The contours correspond to the 20%, 
40%, 60%, and 80% of the maximum line intensity. 


formation of water in the gas phase and is formed via OH'^ 
and H 20 ^ ions, that have recently been found to be impor¬ 
tant tracers of the cosmic-ray ionization rate in the diffuse 
interstellar medium (e.g. Neufeld et al. 2010, Gerin et al. 
2010, Hollenbach et al. 2012). Strong emission lines of OH^, 
H20^, and H 30 ^ ions are also found in AGNs (e.g. Van der 
Werf et al. 2010). Detections of H 30 ^ transitions toward star- 
burst galaxies and AGNs have also been reported with the 
JCMT (e.g. Aalto et al. 2011 and Van der Tak et al. 2008). 

- CS is a well-known tracer of dense gas that has been anal¬ 
ysed both in Galactic star-forming regions (e.g. Helmich & 
van Dishoeck 1997) and in external galaxies, including star- 
burst galaxies (e.g. Aladro et al. 2011b). Multiple rotational 
transitions of CS have previously been studied in W49A by 
Serabyn et al. (1993), leading to the idea that the large star- 
formation activity in W49A has been triggered by a cloud- 
cloud collision. CS J-1-6 has been detected toward W49A 
with a spatial extent of ~ 40" x 24". Its isotopologues '^CS 
and C^"^S have also been detected in the SLS survey. 

- We detected two hydrogen recombination lines: H26a at 
353.6 GHz and H33y6 at 335.2 GHz. Their measured widths 
of ~37 km s 'and ~31 km s 'are larger than the widths 
found for molecular lines detected toward the W49A center 


in the line survey and are similar to those of other H recom¬ 
bination lines detected toward W49A (e.g. Galvan-Madrid 
et al. 2013). The velocity of the H26a' line is consistent 
with that of the H41a line detected by Galvan-Madrid et 
al. (2013). Other hydrogen recombination lines were previ¬ 
ously detected in W49A by Gordon & Walmsley (1990) and 
Galvan-Madrid et al. (2013). 

3.3. Line profiles and velocity structure 

Eigure 8 shows typical line profiles observed toward the vari¬ 
ous subregions of W49A for the most spatially extended molec¬ 
ular lines, with the exception of HCN, HNC, and HCO^ that 
were analysed in Roberts et al. (2011). Most of the detected 
lines are double-peaked or asymmetric. While lines observed 
toward the center and the South-West clump peak around the 
source velocity of ~7-8 km s“' (e.g. Roberts et al. 2011 and 
references therein), line emission is seen mainly blue-shifted to¬ 
ward the Eastern tail, and red-shifted toward the Northern clump, 
similar to what is seen for HNC, HCN, and HCO'^, and their 
isotopologues (Roberts et al. 2011). Some of the asymmetric, 
double-peaked line profiles have a brighter blue-shifted than a 
red-shifted peak, which can be interpreted as gas infalling to- 
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Fig. 8 . Line profiles of molecules observed toward the various subregions of W49A. 


ward the W49A center. Such line profiles that are consistent with 
infall signatures have been observed toward the source center in 
HCO+ 1-0 (Welch et al. 1987, Serabyn et al. 1993) and in the 
SLS survey traced by HCN 4-3 and 3-2, HNC 4-3, and HCO+ 
4-3 (Roberts et al. 201 1). Though most observed lines show an 
excess emission in the blue-shifted wing of the line profile to¬ 
ward the source center, the line profiles which are most consis¬ 


tent with infall signatures are seen in HCN 4-3 and 3-2, HNC 
4-3, and HCO+ 4-3. 

To investigate the velocity distribution of the observed line 
emission, we obtained channel maps for the species shown in 
Fig. 8. CO 3-2 shows emission in a large range of velocities be¬ 
tween -25 and 35 km s ' (Fig. 9). The most extended line emis¬ 
sion is seen around the expected source velocity, in the range 
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Fig. 9. Channel maps of CO at a resolution of '-15" at the given LSR velocities (in km s ^). 
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Fig. 10. Channel maps of '^CO at a resolution of ~15" at the given LSR velocities (in km s '). 


between +5 and +10 km s“'. Toward the Eastern tail, significant 
blue-shifted line emission is seen at velocities between -5 and 
0 km s '. The Northern clump is more prominent at velocities 
of ~2Q km s ', but also shows significant emission around the 
source velocity. The South-west clump is only seen around the 
expected source velocity (between 5 and 10 km s“'). Significant 
line emission is seen at a velocity of +15 km s“' that does not 
correspond to any of the analysed subregions, including the cen¬ 
ter. The '^CO 3-2 line shows a very similar structure (Fig 10) to 
that traced by CO 3-2 with the main difference at a velocity of 
+ 15 km s ', likely related to a self-absorption of the optically 
thick CO 3-2 lines around the center. 


Lines that show spatially extended emission other than CO 
and '^CO are distributed over the velocity range between -5 and 
+20 km s“' (Fig. 11). The strongest emission is seen in CN, 
especially the redshifted wings of the CN lines. While other 
molecules show emission toward the center and the Northern 
clump at a velocity of +15 km s“', CN emission at this velocity 
is also detected toward the South-west clump. 
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Fig. 11. Channel maps of the most spatially extended species in the SLS at the given LSR velocities (in km s '). 
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3.4. Excitation and coiumn densities 


Table 2 includes a number of species with at least three detected 
transitions covering a sufficient energy range to probe excita¬ 
tion conditions. We estimate the excitation conditions of these 
molecules using the population diagram method introduced by 
Goldsmith & Langer (1999). Column densities (A^tot,i) and exci¬ 
tation temperatures (Tex) can be estimated based on 


In 



= ln 


Mou \ 
Grot / 




ln(CG, 


( 1 ) 


where Nu is the observed upper state column density of the 
molecule including line opacity and beam-source coupling ef¬ 
fects; gu is the degeneracy of the upper state; Q^t is the ro¬ 
tational partition function; k is the Boltzmann constant; ft,, is 
the upper level energy; / = ElJQ.a is the source filling fac¬ 
tor with Q„ ~ 15" the beam size and the source size; and 
Ct = t/(1 - e"’’) where r is the optical depth. For a uniform 
beam filling (Q^ ~ Qj) and low optical depth, Eqn. 1 reduces 
to a rotational diagram with a rotational temperature T^ot and to¬ 
tal column density Atot- For a rotational diagram Equation 1 is 
simplified to only the first two terms on the right-hand side of 
the equation. The terms containing the correction for the opacity 
and the beam filling factor are not included. 

We evaluate Eqn. 1 for a set of Ntot.i, Tex, / and Ct (Mot.i, 
Tex)- We apply in the range between 1" and 15" (uniform 
beam filling) and a column density in the range between 10 ^^ 
cm~^ and 10'* cm“^. Eor molecules with just three detected tran¬ 
sitions we use Trot = Tex to derive a best fit source size and col¬ 
umn density. 

The results are summarized in Table 2. The rotation- and 
population diagrams are shown on Eig. 12. The ln{Nu/gu) val¬ 
ues corresponding to the rotational and the population diagrams 
are close to each other, as many of the observed lines are opti¬ 
cally thin (with a few exceptions of some of the SO 2 , CH 3 OH, 
H 2 CO, and **S 02 lines). The scatter in the population diagram 
values is smaller than in the rotational diagrams as the optically 
thicker lines have been corrected for their optical depth. In ad¬ 
dition the lines have been corrected for the best fit size of the 
emitting region. 

Most molecular lines have excitation temperatures between 
-100 and -200 K with H 2 CO, H 2 CS, and SO'*0 showing 
the highest excitation. The population diagrams suggest a non- 
uniform beam filling in several cases and result in emission orig¬ 
inating in regions of -1-4", that are consistent with the sizes of 
hot cores and UCHii regions revealed by high resolution studies 
(e.g. Wilner et al. 2001). Together with the spatially extended 
emission seen in the maps, this is evidence that the gas in W49A 
has a clumpy structure. 

CH 3 CCH has a very low excitation temperature toward the 
center (-32 K) suggesting an origin in colder material around 
the central cluster of UC Hit regions and hot cores. CH 3 CCH 
was also found to originate in a lower excitation component to¬ 
ward Sgr B2, Orion and DR 21 (Churchwell & Hollis 1983). 
The other species that shows a very low excitation temperature 
is (Tex ~ 33 K). Results for SO are not listed in Table 2 as 
the detected SO transitions are not consistent with a single exci¬ 
tation temperature. This may suggest that some of the detected 
SO lines originate in the low excitation component suggested by 
the *^SO lines, and the others from the high-excitation (>100 
K) component that most detected molecular lines originate in. 
Apart from different excitation conditions, gas components with 
different source sizes (beam-filling factors) may contribute to the 
observed line emission, as suggested by the large range of source 


sizes corresponding to the different sulphur-bearing molecules. 
Distinguishing between these explanations requires observations 
of more SO transitions. 

Differences in the excitation of the subregions (Northern 
clump. Eastern tail, and South-west clump) are well traced by 
H 2 CO, as more of its transitions were detected with a spatial 
extent large enough to include the discussed subregions. More 
details are shown in Nagy et al. (2012). 

Eor molecules with no information available on their excita¬ 
tion we derive a column density in the optically thin LTE approx¬ 
imation, assuming rotation temperatures in the range between 75 
K and 150 K, as most molecules with multiple detected transi¬ 
tions have rotation temperatures that fall into this range. Then 
we derive column densities using: 


Ntot 


QjTzod ^EJkT,„, 
hc^ guK\ 



( 2 ) 


Table 3 includes a selection of molecules that have been de¬ 
tected in AGNs, starburst galaxies, and in Galactic star-forming 
regions. The line intensities were obtained using a Gaussian fit 
for nearly Gaussian line profiles. Eor double-peaked or asym¬ 
metric line profiles, the velocity-integrated intensity is used to 
calculate column densities. The error of these column densities 
is about a factor of 2, dominated by the uncertainty in Tex. Eor 
the species where the number of detected transitions was not 
enough to correct for the opacity, the optically thin approxima¬ 
tion results in an under-estimate of the column densities. This is 
the case for some of the species detected toward the center, such 
as for CS and CN. The CS column density is under-estimated by 
a factor of -3, by applying the observed A(C*''S) and an isotopic 
ratio of of 22 (Ererking et al. 1980). The column density 

of CN cannot be constrained as no isotopologue of CN has been 
detected in the SLS survey. Therefore, the value derived from 
CN is a lower limit of the CN column density. We discuss these 
column densities in more details in Sect. 4. 


4. Discussion 

The data reveal a complex structure both in terms of kinematics 
and chemical composition. The molecular line tracers detected 
in the SLS frequency range can be related to various physical 
components including shocked regions and PDRs. Some of the 
detected molecules can be used for a comparison between the 
local starburst seen in W49A and global starburst phenomenon 
seen for some external galaxies. In the next sections we follow¬ 
up on the simple classification of the detected species presented 
in section 3.2 based on kinematics and a comparison to well 
known environments of UV-dominated and shock chemistry. 

4. 1. The kinematics of chemicaiiy reiated species 

In Table A. 1 in Appendix A we list the basic line parameters de¬ 
rived from fitting a single profile Gaussian to each line. Here we 
use the velocities and EWHM line widths to obtain information 
on the kinematics of the previously discussed species. This is a 
simplification, as the fitting of most line profiles requires more 
components than a single Gaussian fit (e.g. Galvan-Madrid et 
al. 2013, Roberts et al. 2011). In this simple comparison fitting 
the line profiles accurately is beyond our scope, which is a com¬ 
parison of the line widths and velocities of the observed lines 
which gives information on the regions where the bulk of emis¬ 
sion arises from. Eor the most asymmetric lines we use only the 
line width derived using a Gaussian fit, and for its vlsr we use 
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Fig. 12. Rotation and population diagram plots toward the center. For CH 3 OH, H 2 CO, and SO 2 the transitions from the high- 
frequency extension of the SLS (between 360-373 GHz) were originally analyzed in Nagy et al. (2012). 


the velocity of the peak intensity. Most of the observed lines cor¬ 
respond to species associated with shock tracers (as explained 
in Sect. 3.2.1) including sulphur-bearing molecules. Fig. 14 in¬ 
cludes a summary of the fitted FWHM and peak velocity values 
observed for the molecules associated with shock tracers. The 
lines of molecules with a possible origin in shocks covers a large 
range of line widths (9-19 km s“*) and peak velocities (7-14 km 
s“'), but most species peak at a velocity of 9-10 km s^'and have 
a width of 15-16 km s '. The large range in the observed pa¬ 
rameters is possibly related to the large region (~ 0.8 pc) covered 
by the JCMT beam in the 345 GHz band, which covers multi¬ 
ple sources, which are only fully resolved with interferometers, 
such as in the maps presented in Galvan-Madrid et al. (2013). 
Figure 13 shows the fitted peak velocities vs the line widths for 
species corresponding to shock tracers (black color), PDR trac¬ 
ers (blue), and complex organic molecules (red). The error bars 
represent the range of the parameters for the different transitions. 
All three groups cover a substantial fraction of the parameter 
space for both velocities and line widths that are discussed in 
this paper. The highest velocities are found for the vibrationally 
excited molecules. For HCN, V 2 =l, the high peak velocities are 
combined with large line widths (FWHM~20 km s '). The fact 
that the molecules which trace different chemistries do not cor¬ 
relate with kinematics may be related to the fact that different 
individual sources are covered by the JCMT beam. The kine- 


matical signatures are likely to be dependent on the source that 
they correspond to, and not on the different chemical groups to 
which they can be related. 


Another method to test whether the assumed chemically re¬ 
lated species indeed trace different regions is a comparison of the 
distributions of the observed line parameters. As most observed 
transitions correspond to shock tracers, we compare the distri¬ 
bution of the observed Vlsr and FWHM line widths to those of 
every other transition. Fig. 14 shows the Vlsr and FWHM line 
width distributions of the species classified as shock tracers and 
of all the other transitions. We compare the distribution of the 
fitted parameters of these groups of species using Kolmogorov- 
Smirnov tests, as shown in Fig. 14. The probability that the two 
distributions are different is about 20 % for the line widths and 
about 38% for the peak velocities. This supports our conclu¬ 
sion that the classification of chemically related species is not 
confirmed by the measured line parameters, likely, because on 
the scale probed in this paper they correspond to many different 
sources, which are not resolved by the JCMT beam. The larger 
difference seen in the velocity distributions compared to the dis¬ 
tributions of the line widths is likely due to the slightly asym¬ 
metric line profiles included in the samples. 
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Table 2. Results of the rotation and population diagram analysis for molecules with multiple transitions. 


Species 

Trot 

Mo, 

Tex 


Source size 

Optical depth^* 


(K) 

(cm“^) 

(K) 

(cm“5) 

(") 


802*’ 

170±23 

(3.1 ±0.02) X 10*** 

118 ±7 

2.1!°;5 X 10*** 

9 Q+0.0 
^•^-0.0 

0 . 1 - 9.3 

33S02 

105±18 

(1.3 ± 0.01) X 10*5 

95!** 

7.6!‘;J X 10*’ 

1.0±0.0 

0 . 8 - 5.3 

3^802 

140±9 

(4.0 ± 0.02) X 10*5 

133±11 

4.8!“ * X 10*5 


0 . 01 - 0.2 

SO 2 , V 2 =l 

146±26 

(1.9 ± 0.06) X 10*’ 

143!“ 

2.0!™ X 10“ 

14 5+^-^ 

0.001-0.004 

CH 3 OH'’ 

75±5 

(7.7 ± 0.04) X 10*5 

65!’ 

5.0!“ X 10*’ 

^•^-0 

0.14-2.69 

H 2 CO'’ 

213±96 

(8.6 ± 0.09) X 10*'* 

159-F90 

1.4!“ X 10*5 

Q "2+4.9 

0.10-1.50 

HNCO 

61 ±20 

(4.9 ± 0.2) X 10"* 

54!S 

7.6!“**** X 10*'* 

12 

0.01-0.05 

H 2 CS 

196±85 

(1.6 ± 0.02) X 10"* 

195!*“ 

1.7!™ X 10"* 

14 

0.001-0.005 

CH 3 CN 

70±30 

(1.1 ±0.1) X 10*5 

69!®^ 

3.1!]5‘* X 10*'* 

- 1.6 

0.16-0.27 

3^80 

31±3 

(5.9 ± 0.04) X 10*'* 

33+8 

6.2![52’ X 10"* 

15.0!“ 

0.01-0.22 

so‘®o 

196±26 

(8.4 ± 0.05) X 10*'* 

1 nn + lKi 

1.91!’;^, X 10“ 

1 q+14.0 

0.20-0.70 

CH 3 CCH 

31±7 

(2.8 ± 0.1) X 10*5 

32 ±7 

2.6!5*4 X 10*5 

15 

i J.U_]2 1 

0.01-0.03 

HC 3 N 

101±39 

(4.7 ± 0.2) X 10*5 

lOU 

4.8!”’ X 10*5 

15 

i^.U_i3 5 

0.004-0.006 

OCS 

92±46 

(1.6 ± 0.08) X 10*5 

92“ 

1.7“53 X 10*5 

1 A 4+0.6 
-13.2 

0.008-0.010 


“ Fixed parameter in the fit 

Transitions in the frequency range between 360 and 373 GHz were included in Nagy et al. (2012). 

The optical depth is calculated for each transition. Here we show the highest and lowest value for each molecule. 

Table 3. Column densities toward the four analysed subregions of W49A (in case of detections) derived in the LTE optically thin 
approximation. The two values correspond to adopted temperatures of 75 K and 150 K, respectively. The error of these column 
densities is about a factor of 2 dominated by the uncertainty in Tex. 


Species 

Transition 

Center 

Column density (cm ’) 

Eastern tail Northern clump 

South-west clump 

C 2 H 

49 / 2 - 37/2 

(1.7-2.8) X 10*5 

(3.8 - 6.0) X 10*+ 

(5.7 - 9.0) X 10*+ 

(5.2 - 8.3) X 10*+ 

NO 

7/2-5/2+ 

(1.0- 1.8) X 10*’ 

so 

0 

X 

1 

00 

0 

(1.0- 1.9) X 10*5 

(4.8 - 8.8) X 10*5 

h*5cn 

4-3 

(1.0- 1.4) X 10*'* 

(6.5 - 9.8) X 10*’ 

(2.0-3.0) X 10“ 

(1.6 - 2.5) X 10“ 

H*5CO+ 

4-3 

(2.5 - 3.8) X 10*5 

(0.8- 1.3) X 10“ 

(1.3 - 2.0) X 10“ 

(1.3 - 1.9) X 10“ 

CN 

M=3-2, y=5/2-3/2 

(1.4-2.2) X 10*5 

(3.0-4.8) X 10*'* 

(4.3 - 6.9) X 10*'* 

(3.7 - 5.9) X 10*+ 

N 2 H+ 

4-3 

(3.0 - 4.4) X 10*5 

(0.7 - 1.1) X 10*5 

(4.1 -6.0) X 10“ 

(1.4-2.0) X 10“ 

cs 

7-6 

(4.4 - 5.6) X 10*'* 

0 

X 

sd 

1 

r-; 

(7.5 - 9.6) X 10*5 

(4.8 - 6.2) X 10*5 

C5'*S 

7-6 

(6.2 - 8.8) X 10*5 

(2.2-2.8) X 10““ 

(4.2 - 6.1) X 10“ 

(5.8-8.3) X 10*' 

H 2 S 

32,1-3i,2 

4.9 X 10*5 




HCO 

4i,4-3i,3 

(1.4-2.4) X 10*5 




CO+ 

yV=3-2, F=5/2-3/2 

(1.4-2.2) X 10*5 




HCS+ 

8-7 

(1.8 - 2.2) X 10*5 




80+ 

8-7 

(4.4 - 5.7) X 10*'* 




H 3 O+ 

32,1-22,0 

(4.8 - 5.3) X 10*'* 




SiO 

8-7 

(0.8- 1.0) X 10"* 





Based on N{CS) and an isotopic ratio of of 22 (Frerking et al. 1980). 


4.2. The importance of UV irradiation 

We probe the importance of UV-irradiation in the chemistry of 
W49A by a comparison of the observed column density ratios 
to values observed toward well-known Galactic PDRs, and by a 
comparison to predictions of PDR models. 

For the comparison of column density ratios to those ob¬ 
served toward PDRs we select a sample of lines that have been 
observed both in the SLS survey and in multiple PDRs (Table 
4). One of the most commonly used tracers of PDRs is the 
[CN]/[HCN] abundance ratio, which was found to be larger than 
unity in NGC 7023 (Fuente et al. 1993). An increase in the abun¬ 
dance ratio of [CN]/[HCN] in PDRs is a result of the production 
of CN via the photodissociation of HCN. Using the HCN col¬ 


umn densities presented in (Nagy et al. 2012) derived in a non- 
FTF approximation and CN column densities from Table 2, we 
estimate a [CN]/[HCN] abundance ratio of ~2.0-4.6 toward the 
Northern clump, ~3.8-6.9 toward the Eastern tail, ~0.4-0.95 to¬ 
ward the South-West clump, and [CN]/[HCN] =0.04-0.07 toward 
the Center. The low [CN]/[HCN] ratio toward the center suggests 
that the chemistry is only partially affected by UV-irradiation. 

Though SiO is a tracer of shocks, it has been detected to¬ 
ward several PDRs as well (e.g. Schilke et al. 2001, Rizzo et al. 
2005). Typical [SiO]/[H'3cO+] abundance ratios of 0.03-0.37 
were found toward Mon R2 (Rizzo et al. 2005), 0.54 toward 
NGC 7023 (Schilke et al. 2001), and 3-14 toward the Orion Bar 
(Schilke et al. 2001), based on detections toward multiple posi- 
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Fig. 14. Top: The distribution of the fitted line widths (left) and peak velocities (right) of the lines corresponding to molecules 
classified as ’shock tracers’. Middle; The distribution of the fitted line widths (left) and peak velocities (right) of every detected line 
in the line survey, except for the double-peaked and very asymmetric line profiles and the shock tracers shown above. Bottom: K-S 
test for the line width (left) and peak velocity (right) distributions shown above, corresponding to shock tracers only and to every 
other molecule that is not classified as a shock tracer. D is the maximum distance between the cumulative fraction of the shock 
tracers and the whole sample of lines and is related to the probability of a difference between the distributions. The part of the 
vertical dotted lines between the two distributions is equivalent to the D parameter. 


tions. We measure an [SiO]/[H'^CO^] abundance ratio of 2.7- 
3.3 toward the W49A center. The HCO abundance with respect 
to the SiO abundance was found to be 50-720 in Mon R2 (Rizzo 
et al. 2005), while is in the range between 16 and 23 toward 
the W49A center. The [HCO]/[SiO] abundance ratio of >4.6- 
10 measured toward NGC 7023 (Schilke et al. 2001) is based 
on an upper limit for the SiO column density. The [HCO]/[SiO] 


abundance ratio toward the Orion Bar was measured to be 2-9, 
based on detections toward multiple positions. Though the large 
(~17 km s ') line width of SiO detected toward the W49A cen¬ 
ter (Table A.l) suggests an origin in shocks rather than PDRs, 
the [HCO]/[SiO] and [SiO]/[H*^CO^] abundance ratios do not 
confirm it. 
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Fig. 13. The measured line widths vs the peak velocities of 
species detected toward the W49A center. The molecules classi¬ 
fied as shock tracers are shown in black, the molecules classified 
as PDR tracers are shown in blue, and the molecules classified 
as complex organic molecules are shown in red. Molecules with 
double-peaked line profiles are not included in this plot. 


The CO^ ion is a tracer of warm (Tgas ^ 500 K) PDR sur¬ 
faces, corresponding to the region of Ay < 2 mag (Fuente et 
al. 2003, Rizzo et al. 2003). We measure an abundance ratio 
of [CO+]/[H'3 cO+] of 0.54-0.56. This is an upper limit on the 
[C0^]/[H'^C0^] ratio as the CO^ line may be blended with 
SO'^O. Ginard et al. (2012) measure [CO+]/[H'3 cO+]= 0.53 in 
Mon R2, similar to that in the Orion Bar (0.52, based on Fuente 
et al. 1996 and Fuente et al. 2003) and to that in NGC 7023 
(0.6, based on Fuente et al. 1993 and Fuente et al. 2003). The 
detection of CO^ at several positions toward the center suggests 
the existence of a hot gas component directly exposed to UV- 
irradiation by a radiation field of Go ~ 3 x 10^ (Vastel et al. 
2001 ). 

HCS^ was found to be abundant in PDRs (Ginard et al. 
2012). We measure a CS abundance with respect to that of 
HCS^ of ~76-108. This is between the lower ratios of 11 and 
25 found toward the ionization front and MP2 positions in Mon 
R2 (Ginard et al. 2012) and the higher [CS]/[HCS+] ratio of 175 
measured toward the Florsehead nebula (Goicoechea et al. 2006). 
This suggests that the [CS]/[HCS^] ratio is not a good tracer 
to probe the contribution of UV irradiation in the chemistry of 
W49A. 

The sulphur-bearing species SO^ and SO 2 can be related to 
shock chemistry, but have also been detected in several PDRs. 
In PDRs SO^ is formed near the PDR surface (at Ay < 2) 
through the reaction - 1 - OH ^ SO^ -t- H and destroyed by dis¬ 
sociative recombination into S and O. SO and SO 2 form in the 
more shielded regions in the PDR. The [SO^]/[S02] abundance 
ratio of <0.02 that we measured toward the center of W49A 
(based on the beam-averaged estimate for SO 2 in Table 2) is at 
least an order of magnitude below the value measured in several 
PDRs, such as in the Orion Bar ([SO^]/[SO2]~0.4-l, Fuente et 
al. 2003). This difference may also be interpreted as the effect 
of shocks on the SO 2 and SO^ the abundances toward the W49 
center. 

The importance of FUV-irradiation in the chemistry of 
W49A can be probed by comparing the observed column den¬ 
sity ratios to predictions of PDR models. Vastel et al. (2001) de¬ 
rive a radiation field toward the W49A center that is equivalent 


Table 4. Summary of observed column density ratios in W49A 
(Center, Eastern Tail, Northern Clump and South-west clump, 
when available) and in Galactic PDRs. References: (a) Fuente 
et al. (1996); (b) Fuente et al. (1993); (c) Ginard et al. (2012); 
(d) Schilke et al. (2001); (e) Rizzo et al. (2005); (f) Fuente et al. 
(2003); (g) Goicoechea et al. (2006) 


Ratio 

W49A 

Orion 

Bar 

NGC 

7023 

Mon 

R2 

Horse- 

head 

CN/HCN center 

0.55 

3.0'“' 

4.5"" 

2-12'“' 


CN/HCN ET 

3.8-6.9 





CN/HCN NC 

2-4.6 





CN/HCN SWC 

0.4-0.95 





SiO/H‘'’CO+ 

2.7-3.3 

3-14<‘‘> 

0.54'‘« 

0.03-0.37<') 


HCO/SiO 

16-23 


>4.6-1 O'"' 

50-720"' 


CO+/H‘’CO+ 

<0.56 

0.52'“'O 

0.6"’"' 

0.53'“' 


CS/HCS+ 

76-108 



11-25"' 

175(8' 

SO+/SO 2 

<0.02 

0.4-1<« 

>0.4(9 

0.2-0.7''' 



to;i^ = 3.5 X 10^;^^o in Draine (1978) units, with;t'o -2.1 x 10“^ 
erg s“* cm“^. We use the 1.4.4 version of the Meudon code (Le 
Petit et al. 2006, Goicoechea &. Le Bourlot 2007, Le Bourlot et 
al. 2012) to calculate column densities for the radiation field de¬ 
rived by Vastel et al. (2001) for expected densities of 10^ cm“^ 
and 10® cm“^ and temperatures of order 100 K derived by pre¬ 
vious studies, corresponding to isobaric models with pressures 
of 10^ cm“^ K and 10^ cm“^ K. We apply a cosmic ray ioniza¬ 
tion rate of 2x10 s ' / H 2 molecule (Hollenbach et al. 2012). 
Predictions of column density ratios of selected species are sum¬ 
marized in Table 5. 

The [HC0]/[H'3C0+] and [C0+]/[H'3C0+] line ratios may 
be consistent with the 10^ cm“^ K pressure PDR model, showing 
that the abundance of some species may be described by PDR 
chemistry. The importance of PDR chemistry toward the region 
in W49A analysed in this paper is further confirmed by Spitzer 8 
fim observations analysed by Peng et al. (2010), showing emis¬ 
sion along arcs and other structures toward most of the region 
where extended emission is seen in PDR tracers such as C 2 H and 
CN. However, the [CN]/[HCN], [CS]/[HCS+], and [SO+]/[S02] 
column density ratios clearly show, that the chemistry of the 
W49A center cannot be explained by a PDR component only, 
and different effects such as shocks also play an important role. 
The orders of magnitude difference seen in the [CS]/[HCS^] and 
[SO'^]/[S 02 ] column density ratios may suggest the importance 
of shock chemistry. The result that the abundances (and abun¬ 
dance ratios) of sulphur-bearing molecules cannot be explained 
by a predominantly PDR chemistry is not surprising. Sulphur 
is released from grains into the gas phase in high temperature 
regions (such as PDRs and hot cores), but can also be evap¬ 
orated by shocks. The regions with different chemistries (and 
different evaporation processes of sulphur from the grains) re¬ 
sult in very different abundances for the sulphur-bearing species. 
This has been observationally confirmed toward regions, where 
the different types of chemistries can more easily be disentan¬ 
gled. One example is Orion KL, where very different abun¬ 
dance ratios of sulphur-bearing molecules have been measured 
toward the Plateau region (dominated by shock chemistry) and 
the Extended ridge (dominated by PDR chemistry) as presented 
in Tercero et al. (2010). 

Based on this comparison we conclude that UV irradiation 
contributes to the chemistry of the W49A center, but it is not the 
dominant effect. 
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Table 5. Comparison of the column density ratios measured 
toward the center of W49A to PDR models with pressures of 
10^ cm“^ K and 10^ cm“^ K illuminated by a radiation field of 
= 3.5 X 10^;ifo- 


Ratio 

W49A center 

Models 

P = 10’ cm-3 K P = 

10* cm-* K 

CN/HCN 

0.55 

42.4 

7.1 

HCO/H'3CO+ 

43.2-75.9 

31.5 

5.0 

CO+/H'3CO+ 

<0.56 

0.81 

0.86 

CS/HCS+ 

76-108 

7151 

5696 

SOTSO 2 

<0.02 

6564 

4892 


4.3. Comparison to regions with shock chemistry 

As we have shown above, the chemistry toward the W49A center 
is mostly related to effects other than FUV-irradiation. Another 
possible effect that has an impact on the chemistry in the W49 
center is shocks related to outflows and winds of the young mas¬ 
sive stars and protostars in the central stellar cluster. The impor¬ 
tance of shocks in W49A has previously been studied by Peng et 
al. (2010) through their relation to expanding shells identified in 
*^CO data, and by Nagy et al. (2012) as a possible contribution 
to the gas heating. Here we select a sample of regions with shock 
chemistry for a comparison to the observed column density ra¬ 
tios (Table 6). 

The Orion KL is a well studied high-mass star-forming re¬ 
gion at a distance of ~414 pc (Menten et al. 2007). Spectral line 
surveys toward Orion KL show a component (Plateau) related 
to outflows and shocks (e.g. Tercero et al. 2010 and references 
therein). We use the OCS/CS, CS/HCS+, CS/HCO+, CS/H 2 CO, 
HCO^/HCS^, and H 2 CO/H 2 CS column density ratios measured 
for the Orion KL plateau by Tercero et al. (2010). The SiO/H 2 CO 
and SO 2 /H 2 CO ratios correspond to the low velocity flow com¬ 
ponent reported by Persson et al. (2007). 

The low-mass Class 0 protostarLl 157 drives a strong molec¬ 
ular outflow (e.g. Bachiller & Perez Gutierrez 1997 and refer¬ 
ences therein). Its blue lobe has a heating rate related to shocks 
comparable to that in the shocked region of Orion KL. The col¬ 
umn density ratios shown in Table 6 correspond to the B1 posi¬ 
tion of the blue lobe measured by Bachiller & Perez Gutierrez 
(1997). 

W3 IRS5 is a bright infra-red source located in the W3 
star-forming region at a distance of 2.0 kpc (e.g. Chavarria et 
al. 2010 and references therein). The embedded young massive 
stars drive several outflows, which affect the chemistry of the re¬ 
gion. In Table 6 we show the column density ratios measured by 
Helmich & van Dishoeck (1997). 

The NGC 1333 star forming region is located at a distance of 
253 pc (Hirota et al. 2008) and contains several low-mass pro¬ 
tostars including IRAS 4A, which is also a source of outflows. 
The column density ratios toward this source in Table 6 are from 
Blake etal. (1995). 

The column density ratios observed toward the W49A cen¬ 
ter are closest to those measured toward the W3 IRS5 source, 
and match within a factor of 3. The lower limit that we derive 
for the CS/HCS^ column density ratio is also close to those to¬ 
ward the Orion KL Plateau and the LI 157 regions. While the 
CS/HCS^ column density is not explained by any of the PDR 
models shown above, they match those observed toward regions 
with shock chemistry reasonably well. The CS/H 2 CO column 
density ratios are also consistent with those measured toward 
LI 157 (within a factor of 2) and NGC 1333 IRAS 4A (within a 
factor of 3). The SiO/H 2 CO and H 2 CO/H 2 CS column density ra- 


Table 6. Comparison of the column density ratios measured 
toward the center of W49A to values measured toward regions 
with shock chemistry. References: Orion KL Plateau SiO/H 2 CO 
and SO 2 /H 2 CO (Persson et al. 2007), Orion KL Plateau other 
column density ratios (Tercero et al. 2010); LI 157 position B1 
(Bachiller & Perez Gutierrez 1997); W3 IRS5 (Helmich & van 
Dishoeck 1997); NGC 1333 IRAS 4A outflow (Blake et al. 
1995). 


Ratio 

W49A 

center 

Orion KL 
Plateau 

L1157 

W3 IRS5 

NGC 1333 
IRAS 4A 

OCS/CS 

0 . 9 - 1.3 

5 

0.2 

0.4 


CS/HCS+ 

76-108 

50 

63 

321 


CS/HCO+ 

0.5-1.0 

5.6 

3.7 

1.1 

19.5 

CS/H 2 CO 

1.6-2.3 

0.1 

0.3-0.9 

2.4 

0.6 

HCO+/HCS+ 

90-165 

9 

17 

304 


H 2 CO/H 2 CS 

5.0-5.1 

18 

2-5.3 

6.2 


SiO/H2CO 

0.1 

0.8 

0.1-0.3 

0.06 

0.04 

S0,/H2C0 

35.4-36.7 

140 

0.4-1 

64.5 

<0.55 


tios are very similar to that measured toward LI 157. The largest 
differences between the regions compared in Table 6 is seen in 
the SO 2 /H 2 CO line ratios. As we have shown in the population 
diagram analysis, the SO 2 lines probably originate in clumps 
smaller than the beam-size (Table 2). Therefore, the large dif¬ 
ferences in the SO 2 /H 2 CO column density ratio between the re¬ 
gions may be related to the combination of the different spa¬ 
tial scales that correspond to the observations and to the clumpy 
structure of SO 2 . 

Even this simple comparison of observed column density ra¬ 
tios shows the importance of shocks in the chemistry of the cen¬ 
ter of W49A. Further evidence could be given by a comparison 
to shock models and by smaller scale observations of sulphur¬ 
bearing molecules in particular. 

4.4. Comparison to starburst gaiaxies and AGNs 

Several molecules seen in our line survey have been detected 
in starburst galaxies and AGNs. A comparison between these 
regions and W49A is shown on Fig. 15, based on line surveys 
toward the starburst galaxy M82 (Aladro et al. 201 lb, carried out 
using the IRAM-30m with beam sizes of 14" - 19" at ~2 mm 
and 9" - 10" at -1.3 mm) and the AGN NGC 1068 (Aladro 
et al. 2013, carried out using the IRAM-30m with a -21" - 29" 
beam). For the column density ratios toward the center the beam- 
averaged column densities are used in the case of species with 
multiple detected transitions, and not the value corresponding to 
the best fit source size indicated by the population diagrams (see 
Sect. 3.4). We use as a reference for the column density 
ratios to be consistent with Aladro et al. (201 lb) and Aladro et 
al. (2013). is also a dense gas tracer and is expected to be 
optically thin. The error bars of the ratios presented in Fig. 15 
are based on the error bars for A(X)/A(C^"^S) for M82 (Aladro 
et al. 201 lb), the column densities for NGC 1068 (Aladro et al. 
2013), and the column densities presented in this paper and in 
Nagy et al. (2012). For some of the abundance ratios this is only 
a lower limit of the error, as those were calculated using column 
density upper limits given for M82 (HNCO and SiO, Aladro et 
al. 2011b) and for NGC 1068 (HjCO, H 2 CS, and SO 2 , Aladro 
et al. 2013). Though a comparison of Galactic scales to those 
probed in external galaxies, such as a scale of -0.8 pc for W49A 
to the 1.5-2 kpc for NGC 1068 (Aladro et al. 2013) and 158-333 
pc for M82 (Aladro et al. 201 lb) is an over-simplification, a few 
conclusions can be drawn. 
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The most similar abundances with respect to that of be¬ 
tween W49A and the AGN and starburst environments are seen 
for C 2 H, HC 3 N, CH 3 CN, CN, H 2 CS, and for H 2 CO in the case 
of M82. In the case of H 2 CO, the Northern clump. Eastern tail, 
and southwest clump regions also show similar abundances w.r.t. 
C^'^S as M82 and the W49A center. C 2 H, CN, and H 2 CO are de¬ 
tected with a large spatial extent in W49A. The similarity of the 
fractional abundance of these species between M82 and W49A, 
and for C 2 H and CN also to that toward NGC 1068 suggests that 
they are good tracers to compare physical and chemical proper¬ 
ties of in Galactic star-forming regions at scales of <1 pc and 
global star-formation seen in external galaxies on scales of >1 
kpc. 

Among the four subregions of W49A, the largest differences 
compared to the other regions in the observed fractional abun¬ 
dances w.r.t. C^"^S are seen for the South-west clump. As C^"^S 
is a dense-gas tracer, the explanation is possibly related to the 
lower average H 2 volume density of the South-west clump com¬ 
pared to the other regions, which is traced by HCN 3-2 and 4-3 
line intensity ratios (Nagy et al. 2012). 

The largest difference compared both to the AGN and star- 
burst examples is seen in the SO 2 fractional abundance. This 
is probably related to the difference between the spatial scales 
of the observations. SO 2 and most sulphur-bearing species in 
W49A are detected in a ~20" x 20" region around the center 
and therefore are least likely to trace ’global’ properties of >kpc 
regions such as the more spatially extended species C 2 H, CN, 
and H 2 CO. Other species with a similar spatial extent to SO 2 
also show significant differences in Fig. 15: CH 3 CCH, HCO, and 
SiO. 

Apart from the abundances, the excitation conditions can 
also be used to compare W49A to external galaxies. We pre¬ 
sented an excitation analysis for 14 molecules, and estimated 
rotational and excitation temperatures. These temperatures can 
also be compared to those found toward M82 and NGC 1068. 
In the following, we use the derived rotational temperatures to 
be consistent with the values derived toward these two galaxies, 
and to use the values that correspond to uniform beam filling 
equivalent to emitting regions of ~0.8 pc (the JCMT beam at 
the observed frequencies). Fig. 16 shows the comparison of the 
species that have a measured excitation temperature toward the 
W49A center and toward M82 (Aladro et al. 201 lb) and toward 
NGC 1068 (Aladro et al. 2013). Most of the rotational tempera¬ 
tures measured toward the W49A center are significantly higher 
(factors of 2-10) compared to those measured toward M82 and 
NGC 1068. This is not surprising due to the very different spatial 
cales that the measurements correspond to. CH 3 CCH is an ex¬ 
ception as the rotational temperature of this molecule very close 
to the value derived toward M82. Apart from CH 3 CCH toward 
M82, CH 3 CN toward M82 and HNCO toward NGC 1068 have 
rotational temperatures close to the value measured toward the 
W49A center. The largest difference is seen for H 2 CO, which 
is not the case when comparing the abundances toward W49A, 
M82, and NGC 1068. This comparison toward W49A is biased 
toward the highest density region with the highest excitation. For 
a better comparison, more transitions toward the off-center re¬ 
gions (such as the Northern clump. Eastern tail, and South-west 
clump) could be compared to similar observations in starburst 
galaxies and AGNs. 

5. Summary 

We presented results from the SFS survey carried out using the 
JCMT at an angular resolution of ~15" in the 330-373 GHz fre¬ 


quency range toward the luminous and massive star-forming re¬ 
gion W49A. Our maps cover a 2x2 arcminutes region around the 
center of W49A, including the regions referred to as the center, 
the Eastern tail, the Northern clump, and the Southwest clump, 
which were selected based on a combination of the morphology 
and the kinematics of the detected molecular lines and were pre¬ 
viously used in Roberts et al. (201 1). 

- The detected 255 transitions correspond to 60 species in¬ 
cluding shock-, and PDR-tracers, and complex organic 
molecules. Excitation conditions can be probed using mul¬ 
tiple detected transitions for 16 species. 

- The species detected with the largest spatial extent include 
CO (and its '^CO and C'^O isotopologues), N 2 H^, and 
H 2 CO; the dense gas tracers including CS, HCN, HNC, and 
HCO^; PDR tracers such as CN and C 2 H. 

- The most spatially extended species reveal a complex spatial 
and kinematic structure, covering the velocity range of [-5, 
h- 20] km s '. Most of the observed lines toward the center are 
double-peaked or asymmetric, with signatures of infall and 
outflow motions. Mostly blue-shifted emission is detected 
toward the Eastern tail region, mostly red-shifted emission 
toward the Northern clump region, while emission toward 
the South-west clump is detected around the source velocity. 

- Based on column density ratios of characteristic species, a 
combination of shock and PDR chemistry affects the physi¬ 
cal and chemical structure of the central 17" x 17" region, 
while UV-irradiation dominates outside of this region includ¬ 
ing the Northern clump. Eastern tail, and South-west clump 
regions. A preliminary comparison to a starburst galaxy and 
to an AGN suggests similar C 2 H, CN, and H 2 CO abundances 
between the ~0.8 pc scale probed for W49A and the >1 kpc 
regions in external galaxies with global star-formation, sug¬ 
gesting them to be best tracers for such a comparison of 
physical and chemical properties of star-formation seen on 
different scales. 

In the future, observations with interferometers, such as with 
ALMA, will probe the chemistry of the W49A region on the 
smaller scales, indicating the differences between the individual 
components of the source that are not resolved by single dish 
telescopes. 
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Fig. 16. A comparison of the rotational temperatures measured toward the W49A center and those measured toward M82 (Aladro 
et al. 201 lb) and toward NGC 1068 (Aladro et al. 2013). 
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Appendix A: The detected lines toward the center of W49A 


Table A.l. The detected species toward the central position of the SLS field. The transitions marked with * have been detected toward AFGL 
2591 in the SLS survey (Van der Wiel 201 1). The transitions marked with * have a non-gaussian line profile. The FWHM line widths of these 
asymmetric lines are from a Gaussian fitting. Their intensity is obtained by integrating over the whole line and their Vlsr is the velocity at Lpeak- 
The parameters without error bars are fixed parameters in the fit. 


Molecule 

Transition 

Freq. 

(MHz) 

£up 

(K) 

(s *) 

J Tmbdv 

(K km s-*) 

Llsr 
( km s-*) 

FWHM 
(km s-*) 

Lpeak 

(K) 

HNCO 

15o,i5-14o,i4 

329664.4 

126.6 

5.04 X 10-^ 

12.09±0.95 

9.68±0.64 

16.57±1.50 

0.69±0.14 

3^802 

82 . 6 -7 1,7 

330191.1 

42.8 

1.19 X 10-“^ 

8.90±0.58 

10.70±0.42 

13.31±1.02 

0.63±0.12 

‘3CO* 

3-2 

330588.0 

31.7 

2.19 X 10-^ 

616.28±1.01 

5.31±0.01 

15.88±0.03 

36.46±1.80 

3^802 

2 l 2 , 20 - 21 l ,21 

330667.6 

218.9 

1.45 X lO-"* 

5.39±0.61 

10.70±0.87 

14.77±1.66 

0.34±0.07 

CHjOH* 

8 - 3 -9-2 E 

330793.9 

146.3 

5.39 X 10-^ 

4.62±0.45 

9.85±0.86 

17.44±1.80 

0.25±0.09 

HNCO 

15i,i4-14i,i3 

330848.6 

170.3 

5.01 X 10-“' 

7.77±0.51 

11.35±0.49 

14.65±1.08 

0.50±0.10 

CHjCN* 

I 83 -I 73 

331014.3 

215.2 

3.07 X 10-3 

2.59±0.25 

9.68±1.00 

10.51±1.00 

0.23±0.11 

CHjCN* 

182-172 

331046.1 

179.5 

3.12 X 10-3 

3.53±0.77 

11.88±1.32 

13.05±2.96 

0.25±0.10 

CHiCN* 

18o-17o 

331071.5 

151.0 

3.16 X 10-3 

9.50±0.26 

11.91±1.00 

17.90±1.00 

0.50±0.11 

CH 3 OH* 

11 ,-lloL^ 

331502.4 

169.0 

3.93 X 10-3 

15.60±0.64 

8.11±0.31 

15.01±0.72 

0.98±0.12 

80. 

1 16,6-125,7 

331580.2 

149.0 

4.35 X 10-3 

17.80±0.61 

10.77±0.24 

14.73±0.61 

1.14±0.12 

80. * 

212 , 20-21 1,21 

332091.4 

219.5 

1.51 X 10-3 

52.14±0.69 

10.08±0.10 

16.28±0.26 

3.01±0.17 

^'^802 

233,21-232,22 

332173.6 

275.1 

2.54 X 10-3 

9.10±0.63 

11.94±0.47 

14.16±1.21 

0.60±0.13 

80. * 

43,1-32,2 

332505.2 

31.3 

3.29 X 10-3 

75.04±0.77 

9.11±0.08 

15.99±0.20 

4.41±0.15 

3^862 

164 , 12 - 163,13 

332836.2 

163.1 

3.02 X 10-3 

17.88±1.02 

12.26±0.44 

16.58±1.23 

1.01±0.13 

8‘"0 

89-78 

333121.6 

75.9 

4.50 X 10-3 

3.89±0.67 

12.38±1.49 

15.91±2.71 

0.23±0.12 

3^80* 

73-67 

333901.0 

79.9 

4.69 X 10-3 

47.01±0.42 

9.57±0.06 

14.00±0.15 

3.15±0.11 

8 O 2 * 

82 , 6 -7 1,7 

334673.4 

43.1 

1.27 X 10-3 

42.19±0.69 

10.38±0.12 

15.60±0.30 

2.54±0.09 

CH 3 OH* 

22-3i a- 

335133.7 

44.7 

2.69 X 10-3 

6.95±0.46 

12.41±0.66 

19.32±1.30 

0.34±0.09 

H 

33/3 

335207.3 



5.19±0.65 

6.62±1.91 

30.50±4.37 

0.16±0.07 

CH 3 CHO 

18o,18,2-17o,17,2 

335318.1 

154.9 

1.30 X 10-3 

6.53±0.41 

11.02±0.61 

20.02±1.44 

0.31±0.06 

CH 3 OH* 

7i-6i A* 

335582.0 

79.0 

1.63 X 10-3 

20.67±0.32 

7.68±0.11 

14.28±0.24 

1.36±0.11 

8 O 2 * 

233,21-232,22 

336089.2 

276.0 

2.67 X 10-3 

37.35±0.41 

10.02±0.08 

14.77±0.20 

2.38±0.10 

HC 3 N 

37-36 

336520.1 

306.9 

3.05 X 10-3 

7.16±0.62 

10.74±0.60 

14.57±1.64 

0.46±0.09 

80 

llio-lOio 

336553.8 

142.9 

6.12 X lO-*^ 

22.03±1.01 

9.40±0.31 

13.57±0.69 

1.52±0.08 

80. 

167 , 9 - 175,12 

336669.6 

245.1 

5.84 X 10-3 

16.22±0.11 

11.05±0.17 

11.46±0.34 

1.33±0.14 

CH 3 OH* 

12 i- 12 o A=" 

336865.1 

197.1 

4.07 X 10-3 

14.71 ±0.44 

7.69±0.22 

15.21±0.53 

0.91±0.10 

c‘"o* 

3-2 

337062.0 

32.4 

4.30 X 10-3 

84.53±0.39 

7.39±0.03 

13.62±0.07 

5.83±0.32 

3380 

73-67 

337199.4 

80.5 

4.65 X 10-3 

19.13±0.64 

10.11 ±0.24 

14.68±0.59 

1 . 22 ± 0.10 

HC3N,V7 = 1 

7=37-36 

337344.7 

628.5 

3.05 X 10-3 

3.19±0.55 

14.02±1.04 

11.72±2.43 

0.26±0.09 

C3^8* 

7-6 

337396.5 

50.2 

8.00 X 10-3 

24.98±0.56 

8.02±0.15 

13.79±0.35 

1.70±0.15 

3^80* 

83-77 

337580.1 

86.1 

4.89 X 10-3 

46.74±0.43 

9.08±0.06 

13.66±0.14 

3.21±0.10 

CH3OH, v,=l 

77,1 - 66,1 A* 

337643.9 

365.4 

1.69 X 10-3 

3.54±0.50 

9.49±1.20 

15.61±2.04 

0.21±0.06 

CH 30 CH 3 

74 , 3,0 - 63 , 4,0 

337787.2 

48.0 

1.94 X 10-3 

1.80±0.53 

11.72±1.38 

10.13±3.81 

0.17±0.04 

3^80 

33-23 

337892.2 

25.3 

1.40 X 10-3 

4.59±0.54 

10.81±0.70 

12.32±1.81 

0.35±0.08 

H2C8 

10i,io-9i,9 

338083.2 

102.4 

5.77 X 10-3 

7.31±0.57 

7.36±0.54 

13.65±1.17 

0.50±0.07 

CH 3 OH* 

7o-6o E 

338124.5 

78.1 

1.70 X 10-3 

15.87±0.88 

6.64±0.37 

13.31±0.85 

1 . 12 ± 0.10 

8 O 2 * 

184 , 14 - 183,15 

338306 

196.8 

3.27 X 10-3 

29.82±0.39 

7.47±0.11 

18.52±0.28 

1.51±0.10 

CH 3 OH* 

7_i-6_i E 

338344.6 

70.6 

1.67 X 10-3 

23.99±0.31 

6.59±0.09 

13.53±0.21 

1.67±0.10 

CH 3 OH* 

7o-6o A* 

338408.7 

65.0 

1.70 X 10-3 

22.13±0.30 

6.62±0.08 

12.26±0.19 

1.70±0.10 

CH 3 OH* 

74,4 — 64,3 A 

338512.6 

145.3 

1.15 X 10-3 

10.32±0.48 

7.35±0.32 

13.83±0.72 

0.70±0.09 

CH 3 OH 

73,4 “ 63,3 A 

338543.2 

114.8 

1.39 X 10-3 

10.68±0.50 

7.97±0.33 

13.87±0.75 

0.72±0.09 

80. * 

2O1,19-192,18 

338611.8 

198.9 

2.87 X 10-3 

40.93±0.50 

6.95±0.09 

15.73±0.22 

2.44±0.10 

CH 3 OH 

72 - 62 A+ 

338639.9 

102.7 

1.57 X 10-3 

7.93±0.90 

6.75±1.00 

11.79±1.00 

0.63±0.10 

CH 3 OH* 

72-62 E 

338721.6 

87.3 

1.55 X 10-3 

21.52±0.27 

5.92±0.09 

13.16±0.20 

1.54±0.17 

3^80. 

144 , 10 - 143,11 

338785.7 

134.5 

3.08 X 10-3 

5.35±0.32 

11.23±0.38 

12.57±0.83 

0.40±0.04 

80* 

33-32 

339341.5 

25.5 

1.45 X 10-3 

18.61±0.34 

8 . 02 ± 0.12 

13.74±0.28 

1.27±0.08 

33802 

164 , 12 - 163,13 

339482.3 

166.1 

3.18 X 10-3 

2.72±0.31 

10.53±0.71 

12.33±1.53 

0.21±0.05 

CN* 

A=3-2, 7=5/2-5/2 

339516.6 

32.6 

2.54 X 10-3 

3.15±0.50 

7.46±1.40 

17.27±3.49 

0.17±0.07 

3'‘80* 

93-87 

339857.3 

77.3 

5.08 X 10-3 

70.69±0.33 

9.15±0.03 

14.43±0.08 

4.60±0.15 

CN** 

A=3-2, 7=5/2-3/2 

340035.4 

32.6 

2.89 X 10-3 

51.50±1.01 

6.98±1.0 

17.38±0.68 

3.07±0.06 

CH 3 OH* 

22-3i a* 

340141.2 

44.7 

2.78 X 10-3 

4.38±0.35 

8.03±0.60 

14.91±1.39 

0.28±0.05 

CN** 

A=3-2, 7=7/2-5/2 

340248.5 

32.7 

3.67 X 10-3 

51.14±0.99 

5.44±1.0 

13.25±1.42 

3.61±0.07 

‘ 3 CH 3 OH 

22,1 - 3i,2 a 

340313.9 

44.6 

2.82 X 10-3 

31.79±0.32 

8.21±0.07 

14.43±0.17 

2.07±0.08 

OC 8 * 

28-27 

340449.3 

237.0 

1.15 X 10-3 

10.68±0.41 

10.38±0.29 

14.55±0.58 

0.69±0.08 

33802 

2O1,19-192,18 

340526.1 

201.5 

2.96 X 10-3 

4.61±0.32 

11.27±0.40 

11.24±0.86 

0.39±0.06 

HC‘30+ 

4-3 

340630.7 

40.9 

3.11 X 10-3 

6.98±0.36 

7.58±0.29 

11.26±0.68 

0.58±0.07 


Continued on next page 
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Table A.l. The detected species toward the central position of the SLS field. The transitions marked with * have been detected toward AFGL 
2591 in the SLS survey (Van der Wiel 2011). The transitions marked with * have a non-gaussian line profile. The FWHM line widths of these 
asymmetric lines are from a Gaussian fitting. Their intensity is obtained by integrating over the whole line and their Vlsr is the velocity at Tpeak- 
The parameters without error bars are fixed parameters in the fit. 


Molecule 

Transition 

Freq. 

(MHz) 

Fup 

(K) 

(s ') 



J Finbdv 

(K km s-*) 

Llsr 
( km s-') 

FWHM 
(km s-*) 

Lpeak 

(K) 

SO* 

78-67 

340714.2 

81.2 

4.99 

X 

Kf^ 

213.26±1.09 

7.49±0.04 

15.61±0.10 

12.83±0.52 

«so 

88-77 

340839.6 

86.8 

5.03 

X 

lO-"* 

22.64±0.26 

10.49±0.08 

13.84±0.18 

1.54±0.09 

so. 

218.14 - 227,15 

341275.5 

369.1 

6.86 

X 

10-5 

14.77±0.36 

11.60±0.15 

12.69±0.37 

1.09±0.10 

HCS+ 

8-7 

341350.2 

73.7 

8.35 

X 

10-* 

5.68±0.37 

9.40±0.46 

13.38±0.87 

0.40±0.08 

CHiOH* 

7i-6i A- 

341415.6 

80.1 

1.71 

X 

lO-'^ 

36.71±0.48 

10.62±0.14 

20.43±0.31 

1.69±0.19 

HCO 

4i,4 - 3],3 

341671.6 

73.8 

3.46 

X 

10-* 

18.60±0.31 

7.78±0.14 

16.70±0.33 

1.05±0.09 

33SO, 

13,,12 - 12i,ii 

341721.7 

94.1 

2.30 

X 

lO-'^ 

10.70±1.45 

10.51±0.95 

17.67±2.28 

0.57±0.09 

CHiCCH 

20-19 

341741.0 

172.2 

1.39 

X 

lO-'^ 

9.38±1.43 

9.82±0.86 

14.53±1.58 

0.61±0.09 

^'^SOz 

53,3 - 42,2 

342208.9 

35.1 

3.10 

X 

10-^ 

15.29±0.51 

11.77±0.19 

11.83±0.46 

1.21±0.10 

3^SO, 

20],19 - 192,18 

342231.6 

198.5 

3.06 

X 

10-* 

16.20±0.52 

11.60±0.21 

14.02±0.55 

1.09±0.10 

^^SOj 

124,8 - 123,9 

342332.0 

109.7 

3.06 

X 

10-* 

15.07±0.47 

11.57±0.20 

12.77±0.46 

1.11±0.10 

CH30H* 

13i-13o A* 

342729.8 

227.5 

4.23 

X 

10-* 

9.95±0.38 

9.06±0.26 

13.71±0.57 

0.68±0.09 

O'^C^'^S 

29-28 

342759.8 

246.8 

1.18 

X 

10-* 

17.47±0.35 

9.66±0.14 

13.45±0.34 

1.22±0.09 

cs* 

7-6 

342882.9 

65.8 

8.40 

X 

10-* 

155.79±0.66 

6.35±0.03 

14.21±0.08 

10.30±0.80 

H2CS 

lOo.io - 9o,9 

342946.4 

90.6 

6.08 

X 

10-* 

4.15±0.35 

9.11±0.52 

11.58±0.99 

0.34±0.06 

HO'^C* 

4-3 

342983.3 

41.2 

3.34 

X 

10-5 

2.49±0.35 

10.83±0.82 

11.91±2.29 

0.20±0.07 

^^SO 

98 - 87, F = 

343087.3 

78.0 

5.09 

X 

10-* 

30.65±0.53 

9.66±0.12 

14.21±0.29 

2.03±0.10 

HfCO 

5i5-4i4 

343325.7 

61.3 

1.12 

X 

10-5 

8.76±0.58 

10.88±0.60 

18.27±1.50 

0.45±0.07 

H2CS 


343414.1 

209.1 

5.56 

X 

10-* 

5.62±0.33 

12.12±0.48 

15.75±0.89 

0.33±0.09 

H2CS 

102,8 ~ 92,7 

343813.2 

143.4 

5.88 

X 

10-* 

2.56±0.41 

10.25±0.83 

11.91±2.68 

0.20±0.10 

HC'^N* 

4-3 

344200.1 

41.3 

1.88 

X 

10-5 

13.85±0.46 

8.49±0.24 

14.51±0.56 

0.90±0.09 

3^S02 

104,6 —103,7 

344245.3 

88.5 

2.96 

X 

10-* 

12.41±0.36 

10.85±0.17 

11.56±0.39 

1.01±0.07 

SO* 

88-77 

344310.6 

87.5 

5.19 

X 

10-* 

199.46±1.96 

7.39±0.07 

15.37±0.17 

12.19±0.50 

^'‘SO, * 

19i,i9-18o,i8 

344581.0 

167.7 

5.16 

X 

10-* 

20.86±0.43 

11.35±0.12 

12.20±0.30 

1.61±0.15 

^^S02 * 

134,10-133,11 

344807.9 

121.6 

3.17 

X 

10-* 

14.45±0.54 

11.54±0.25 

13.69±0.64 

0.99±0.12 

3'‘S02 

154,12-153,13 

344987.6 

148.3 

3.27 

X 

10-* 

26.22±0.64 

6.18±0.22 

17.92±0.49 

1.37±0.14 

33SO, 

144,10-43,11 

345134.4 

137.1 

3.19 

X 

10-® 

7.92±1.14 

12.86±0.93 

15.36±2.61 

0.48±0.12 

so. 

55,1-64,2 

345149.0 

75.1 

9.81 

X 

lo-** 

8.51±1.20 

11.36±0.41 

9.70±1.13 

0.82±0.12 

3^s6, 

84,4-83,5 

345168.7 

71.0 

2.75 

X 

10-* 

12.42±0.68 

12.25±0.28 

11.39±0.83 

1.02±0.12 

^^S02 

94,6-93,7 

345285.6 

79.3 

2.88 

X 

10-* 

13.67±0.66 

11.93±0.25 

10.22±0.56 

1.26±0.17 

H'^CN* 

4-3 

345339.8 

41.4 

1.90 

X 

10-5 

110.90±0.72 

8.59±0.05 

15.70±0.12 

6.64±0.46 

so. 

269,17 - 278,20 

345449.0 

521.0 

7.63 

X 

10-5 

11.82±0.65 

11.12±0.35 

13.25±0.90 

0.84±0.17 

3'‘s62 * 

74,4-73,5 

345519.7 

63.7 

2.59 

X 

10-* 

14.84±0.63 

11.70±0.24 

11.80±0.62 

1.18±0.13 

^'*S02 

64,2-63,3 

345553.1 

57.3 

2.35 

X 

10-* 

11.10±2.13 

12.08±0.96 

10.04±2.15 

1.04±0.13 

^3S02 

19i,i9 - 18o,i8 

345584.7 

170.2 

5.19 

X 

10-* 

4.22±0.60 

12.85±0.53 

7.71±1.28 

0.51±0.11 

HC3N 

38-37 

345609.0 

323.5 

3.30 

X 

10-5 

6.83±0.41 

12.06±0.35 

13.19±1.07 

0.49±0.09 

^'*S02 

54,2-53,3 

345651.3 

51.8 

1.96 

X 

10-* 

6.12±0.78 

11.99±0.51 

7.46±0.98 

0.77±0.11 

^^S02 

44,0 - 43,1 

345678.8 

47.2 

1.31 

X 

10-* 

3.38±0.39 

11.30±0.61 

10.38±1.30 

0.31±0.09 

CO* 

3-2 

345796.0 

33.2 

2.50 

X 

lo-** 

1161.6±4.43 

2.77±0.001 

21.22±0.10 

51.42±7.3 

3^SO, 

174,14-173,15 

345929.3 

178.8 

3.37 

X 

10-* 

9.96±0.53 

12.02±0.32 

12.24±0.77 

0.76±0.09 

NS* 

15/2-13/2 

346220.1 

71.0 

7.38 

X 

10-* 

7.35±0.50 

10.84±0.68 

20.31 ±1.54 

0.34±0.06 

SO2, V2=l 

19i,i9 - 18o,i8 

346379.2 

930.6 

5.16 

X 

10-* 

9.82±0.41 

9.03±0.38 

18.64±0.91 

0.49±0.09 

HC3N, V7=l 

7=38-37 

346455.7 

645.1 

3.31 

X 

10-5 

2.10±0.34 

13.99±0.96 

9.77±2.07 

0.20±0.07 

so* 

98-87 

346528.5 

78.8 

5.38 

X 

10-* 

231.50±0.69 

7.63±0.02 

16.69±0.06 

13.03±0.68 

”S02 

53,3-42,2 

346589.8 

36.0 

1.76 

X 

10-5 

6.41±0.33 

10.91±0.30 

11.76±0.71 

0.51±0.06 

S02 * 

19i,i9-18o,i8 

346652.2 

168.1 

5.22 

X 

10-* 

96.43±0.32 

9.42±0.03 

15.47±0.06 

5.85±0.14 

HC3N V7=l 

7=38-37 

346949.1 

645.6 

3.32 

X 

10-5 

1.55±0.30 

13.93±0.73 

6.88±1.66 

0.21±0.09 

H'^co* * 

4-3 

346998.3 

41.6 

3.29 

X 

10-5 

50.66±0.41 

6.58±0.05 

11.86±0.12 

4.01±0.19 

SiO 

8-7 

347330.6 

75.0 

2.20 

X 

10-5 

68.67±0.66 

7.91±0.08 

16.98±0.20 

3.80±0.16 

3'*S02 

282,26-281,27 

347483.1 

391.2 

2.65 

X 

10-* 

5.01±0.37 

11.48±0.45 

11.89±1.03 

0.40±0.07 

so* 

8-7 

347740 

70.1 

2.28 

X 

10-* 

19.44±0.41 

8.83±0.15 

13.87±0.33 

1.32±0.09 

SO,,v,=l 

132,12 - 12l,ll 

347991.8 

854.3 

2.41 

X 

10-* 

4.48±0.47 

9.98±0.88 

16.43±1.81 

0.26±0.08 

3'*S02 

194,16-193,17 

348117.5 

212.9 

3.50 

X 

10-* 

28.62±0.38 

11.22±0.09 

14.24±0.22 

1.89±0.10 

HC'^O* 

4-3 

348211.2 

41.8 

3.32 

X 

10-5 

1.64±0.33 

8.75±1.12 

10.73±2.59 

0.14±0.05 

HN'^C 

4-3 

348340.8 

41.8 

2.03 

X 

10-5 

5.66±1.05 

11.00±1.00 

11.88±1.00 

0.45±0.11 

SO, * 

242 , 22 - 233,21 

348387.8 

292.7 

1.91 

X 

10-* 

48.45±1.14 

10.12±1.00 

13.99±1.00 

3.25±0.10 

33s6, 

124,8-123,9 

348490.9 

111.9 

2.96 

X 

10-® 

6.28±0.78 

8.90±1.17 

18.08±3.31 

0.33±0.10 

0-H2CS* 

10i,9-9i,8 

348534.4 

105.2 

6.32 

X 

10-* 

13.13±0.62 

9.32±0.36 

15.55±0.86 

0.79±0.10 

C2H* 

49/2-35/2, F=4-3 

349108.5 

41.9 

5.76 

X 

10-* 

9.05±0.46 

9.59±0.37 

14.35±0.80 

0.59±0.08 

q 

00 

0 

154,12-153,13 

349209.6 

149.8 

3.37 

X 

10-* 

2.84±0.62 

10.65±1.58 

12.26±3.04 

0.22±0.09 


Continued on next page 
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Table A.l. The detected species toward the central position of the SLS field. The transitions marked with * have been detected toward AFGL 
2591 in the SLS survey (Van der Wiel 2011). The transitions marked with * have a non-gaussian line profile. The FWHM line widths of these 
asymmetric lines are from a Gaussian fitting. Their intensity is obtained by integrating over the whole line and their Vlsr is the velocity at Tpeak- 
The parameters without error bars are fixed parameters in the fit. 


Molecule 

Transition 

Freq. 

(MHz) 

£up 

(K) 

(s *) 



J L,„bdv 

(K km s-*) 

Llsr 
( km s-') 

FWHM 
(km s-*) 

Lpeak 

(K) 

S0‘“0 

134,10-133,11 

349224.9 

123.3 

3.27 

X 

Kf^ 

2.44±0.53 

11.95±0.28 

4.76±1.27 

0.48±0.09 

CzH* 

49 / 2 - 37/2 

349338.3 

41.9 

1.28 

X 

10-* 

33.15±0.46 

8.49±0.10 

13.97±0.22 

2.23±0.17 

CzH* 

47 / 2 - 35/2 

349400.7 

41.9 

1.20 

X 

10-^ 

28.34±1.14 

10.32±0.24 

14.04±0.46 

1.90±0.10 

CHjCN* 

192-182 

349426.8 

196.3 

3.68 

X 

10-3 

5.91±3.13 

10.57±2.57 

19.45±12.17 

0.29±0.10 

CHiCN* 

19o-18o 

349453.7 

167.7 

3.72 

X 

10-3 

9.74±1.63 

12.05±1.12 

15.03±1.74 

0.61±0.10 

”S02 

104,6-103,7 

350302.4 

90.4 

4.02 

X 

10-*^ 

2.87±0.34 

10.39±0.73 

11.41±1.29 

0.24±0.08 

HNCO* 

16i,16-15i,15 

350333.1 

186.2 

5.97 

X 

10-* 

5.82±0.32 

12.30±0.26 

9.89±0.64 

0.55±0.08 

NO* 

7/2-5/2+ 

350694.8 

36.1 

4.81 

X 

10-'^ 

35.27±0.44 

12.03±0.09 

14.78±0.21 

2.24±0.06 

33SO, 

134 , 10 - 133,11 

350787.0 

124.0 

2.64 

X 

10-® 

2.84±0.36 

11.42±0.89 

14.02±2.08 

0.19±0.08 

so, * 

106,4-1 15,7 

350862.8 

138.8 

4.40 

X 

10-3 

24.44±0.51 

11.46±0.13 

12.58±0.30 

1.83±0.12 

CHjOH* 

li-OoA* 

350905.1 

16.8 

3.31 

X 

10-^ 

28.53±0.52 

5.82±0.11 

12.22±0.28 

2.19±0.24 

SO'^O 

124 , 8 - 123,9 

350993.2 

112.1 

3.27 

X 

10-3 

5.55±0.54 

11.15±0.46 

9.97±1.17 

0.52±0.12 

NO* 

7/2-5/2- 

351049.0 

36.1 

4.99 

X 

10-3 

35.03±0.61 

9.79±0.14 

16.80±0.34 

1.96±0.13 

^^SOj 

84 , 4 - 83,5 

351177.4 

72.7 

5.65 

X 

10-3 

4.14±0.43 

12.48±0.38 

8.01±1.05 

0.49±0.09 

so, * 

53,3-42,2 

351257.2 

35.9 

3.36 

X 

10-3 

42.96±0.59 

9.27±0.11 

16.27±0.28 

2.48±0.13 

33s6, 

24 , 4 - 73,5 

351508.5 

65.2 

6.76 

X 

10-3 

3.09±0.34 

11.35±0.47 

8.20±0.96 

0.35±0.09 

”S 02 

64 , 2 - 63,3 

351542.1 

58.7 

8.06 

X 

10-3 

8.37±0.55 

13.22±0.61 

18.63±1.30 

0.42±0.09 

HNCO 

16o,i6 “ 15o,i5 

351633.3 

143.5 

6.13 

X 

10-3 

20.41±0.79 

12.36±0.24 

12.78±0.57 

1.50±0.12 

”S 02 

44 , 0 - 43,1 

351661.8 

48.5 

1.18 

X 

10-3 

3.39±0.53 

9.73±0.85 

10.67±1.90 

0.30±0.11 

33SO, 

174 , 14 - 173,15 

351743.7 

181.9 

1.78 

X 

10-3 

5.17±0.53 

13.39±0.56 

11.22±1.08 

0.43±0.18 

0 -H 2 CO* 

5i,5-4i,4 

351768.6 

62.5 

1.20 

X 

10-3 

104.71 ±1.35 

6.54±0.08 

13.75±0.22 

7.15±0.65 

so, * 

144 , 10 - 143,11 

351873.9 

135.9 

3.43 

X 

10-3 

63.49±0.71 

10.17±0.08 

14.94±0.19 

3.99±0.17 

3^s6, 

214,18-213,19 

352082.9 

250.8 

3.65 

X 

10-3 

8.53±0.59 

11.65±0.37 

11.33±0.96 

0.71 ±0.07 

ocs* 

29-28 

352599.6 

253.9 

1.28 

X 

10-3 

11.23±0.33 

10.91±0.20 

13.54±0.44 

0.78±0.08 

HNCO* 

16i,i5-15i,i4 

352897.6 

187.2 

6.10 

X 

10-3 

6.28±0.33 

11.71 ±0.27 

10.12±0.63 

0.58±0.07 

3^SO, 

147,7-156,10 

353002.4 

212.6 

5.70 

X 

10-3 

1.18±0.20 

11.24±0.98 

10.59±1.51 

0.11 ±0.05 

so‘*o 

263 , 24 - 262,25 

353195.8 

342.7 

2.94 

X 

10-3 

1.59±0.25 

12.54±0.89 

11.07±1.87 

0.13±0.03 

H 

26a 

353622.8 





18.97±0.46 

11.63±0.44 

37.31±1.07 

0.48±0.06 

CO*^ 

N=3-2, F=f-f 

353741.3 

33.9 

2.06 

X 

10-3 

5.55±0.45 

12.11±0.82 

21.09±2.06 

0.25±0.06 

HCN,V 2 = 1 ““‘‘/' * 

4-3 

354460.4 

1066.9 

1.87 

X 

10-3 

9.78±0.95 

12.50±1.0 

15.00±1.00 

0.62±0.06 

HCN*’*'^"'*^^ *♦ 

4-3 

354505.5 

42.5 

2.05 

X 

10-3 

184.9±1.30 

1.32±1.0 

22 . 00 ± 2.00 

8.55±0.06 

SO 2 , V 2 =l 

164,12 - 163,13 

354800 

927.8 

3.60 

X 

10-3 

3.16±0.32 

10.73±0.91 

15.39±1.89 

0.19±0.08 

SO, * 

124,8 - 123,9 

355045.5 

111.0 

3.40 

X 

10-3 

49.46±0.32 

9.74±0.05 

15.04±0.12 

3.09±0.09 

S 02 

174,14 - 18i,i7 

355186.5 

180.1 

2.62 

X 

10-3 

5.16±0.36 

9.08±0.57 

16.41±1.47 

0.30±0.05 

H'^NC 

4-3 

355439.5 

42.6 

1.69 

X 

10-3 

1.25±0.23 

10.68±0.83 

8.36±1.71 

0.14±0.05 

s‘®o 

89 - Is 

355571.1 

93.1 

5.74 

X 

10-3 

8.35±0.26 

11.21±0.24 

16.00±0.58 

0.49±0.06 

CH 2 OH* 

13o,i3-12i,i2 a* 

355603.1 

211.0 

2.53 

X 

10-3 

9.52±0.27 

9.61±0.22 

15.04±0.47 

0.59±0.06 

CH 2 OH* 

15i,i4-15o,i5 a* 

356007.2 

295.3 

4.60 

X 

10-3 

7.87±0.62 

8.70±0.58 

14.19±1.13 

0.52±0.10 

SO, 

157,9 “ 166,10 

356040.6 

230.4 

6.40 

X 

10-3 

23.58±0.67 

10.51±0.19 

13.27±0.45 

1.67±0.10 

3^s62 

253,23 “ 252,24 

356222.2 

320.0 

2.97 

X 

10-3 

6.39±0.51 

10.72±0.45 

10.87±0.92 

0.55±0.08 

HCN, V 2 = 1 

7 = 4-3 

356255.6 

1067.12 

5.33 

X 

10-3 

11.42±0.71 

12.20±0.57 

20.00±1.61 

0.54±0.07 

HCO+ * 

4-3 

356734.2 

42.8 

3.57 

X 

10-3 

503.55±0.55 

3.10±0.01 

27.35±0.03 

17.30±1.30 

3^SO, * 

20o,2o - 19i,i9 

357102.2 

184.8 

5.81 

X 

10-3 

25.19±0.53 

11.25±0.14 

13.36±0.33 

1.77±0.14 

so, * 

134,10 “ 133,11 

357165.4 

123.0 

3.51 

X 

10-3 

72.28±0.41 

9.24±0.04 

15.20±0.10 

4.47±0.15 

so, * 

154,12 - 153,13 

357241.2 

149.7 

3.62 

X 

10-3 

67.49±0.68 

9.24±0.07 

14.70±0.17 

4.31±0.14 

so, * 

1 14,8 - 113,9 

357387.6 

100.0 

3.38 

X 

10-3 

71.00±0.84 

9.85±0.09 

15.00±0.21 

4.45±0.13 

3'‘s62 

325,27 - 324,28 

357497.8 

547.5 

4.55 

X 

10-3 

2.34±0.24 

12.21±0.47 

9.23±1.22 

0.24±0.05 

S 02 * 

84,4 - 83,5 

357581.4 

72.4 

3.06 

X 

10-3 

81.60±0.61 

9.32±0.06 

17.96±0.16 

4.27±0.24 

S 02 * 

94,6 - 93,7 

357671.8 

80.6 

3.20 

X 

10-3 

86.02±0.30 

9.88±0.03 

16.93±0.07 

4.77±0.23 

go,!)!™*)/! * 

74,4 - 73,5 

357892.4 

65.0 

2.87 

X 

10-3 

78.52±2.36 

9.36±1.00 

15.65±1.00 

4.71±0.15 

g0_^blend/2 ★ 

64.2 “ 63,3 

357925.8 

58.6 

2.60 

X 

10-3 

74.40±0.36 

9.47±0.04 

15.64±0.09 

4.47±0.14 

gO,blcnii/3 

174,14 - 173,15 

357962.9 

180.1 

3.73 

X 

10-3 

69.61±0.38 

10.06±0.04 

15.86±0.10 

4.12±0.14 

g0_^bleiid/4 ★ 

54,2 - 53,3 

358013.2 

53.1 

2.18 

X 

10-3 

68.11±0.36 

9.71±0.01 

15.47±0.10 

4.14±0.14 

gO,blcn‘i/5 ♦ 

44,0 - 43,1 

358037.9 

48.5 

1.45 

X 

10-3 

48.82±2.39 

9.68±1.00 

13.32±1.00 

3.44±0.14 

SO, * 

20o,2o - 19i,i9 

358215.6 

185.3 

5.83 

X 

10-3 

102.21±0.56 

9.29±0.04 

15.47±0.10 

6.21±0.17 

^^s 62 

234,20 - 233,21 

358347.3 

292.4 

3.86 

X 

10-3 

9.32±0.30 

11.67±0.20 

12.70±0.48 

0.69±0.08 

CH 3 OCH 3 

55,1,1 - 44,1,1 

358454.0 

48.8 

2.94 

X 

10-3 

3.32±0.30 

8.56±0.68 

14.95±1.40 

0.21±0.05 

CH 3 OH* 

4i - 3o £ 

358605.8 

44.3 

1.32 

X 

10-3 

24.58±1.35 

7.66±0.38 

13.84±0.83 

1.67±0.14 


Continued on next page 
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Table A.l. The detected species toward the central position of the SLS field. The transitions marked with * have been detected toward AFGL 
2591 in the SLS survey (Van der Wiel 201 1). The transitions marked with * have a non-gaussian line profile. The FWHM line widths of these 
asymmetric lines are from a Gaussian fitting. Their intensity is obtained by integrating over the whole line and their Vlsr is the velocity at Tpeak- 
The parameters without error bars are fixed parameters in the fit. 


Molecule 

Transition 

Freq. 

(MHz) 

£up 

(K) 

(s ') 



J Lmbdv 

(K km s-*) 

Llsr 
( km s-') 

FWHM 
(km s-*) 

Lpeak 

(K) 


% - 85 

358645.7 

99.3 

5.92 

X 

Kf^ 

9.23±3.28 

9.38±2.27 

12.81±5.15 

0.68±0.05 


2 I 3 - 2 O 3 

358756.5 

254.5 

1.58 

X 

lO-"* 

3.63±0.18 

9.34±0.30 

12.03±0.66 

0.28±0.04 

CH3CCH'’'“‘*^ 

2 I 2 - 2 O 2 

358790.6 

218.3 

1.60 

X 

10 -^ 

2.75±0.28 

8.37±0.67 

12.40±1.45 

0.21±0.06 

CH3CCH'’'“‘‘/3 

2 I 0 - 2 O 0 

358817.9 

189.5 

1.61 

X 

lO-"* 

7.92±0.43 

10.53±0.36 

14.48±0.93 

0.51±0.06 

so‘’o^ 

152,14-14i,13 

358986.2 

119.9 

2.84 

X 

10 -^ 

20.64±0.25 

9.95±0.08 

13.18±0.19 

1.47±0.10 

S 02 * 

253,23 - 252,24 

359151.2 

320.9 

3.10 

X 

10 -“* 

42.63±0.21 

10.57±0.04 

14.47±0.09 

2.77±0.08 

3^802 

242,22 “ 233,21 

359651.7 

292.4 

2.20 

X 

10 -“ 

6.76±0.34 

12.02±0.26 

11.01±0.71 

0.58±0.07 

S 02 * 

194,16 - 193,17 

359770.7 

214.3 

3.85 

X 

10 -“ 

58.89±0.28 

10.13±0.03 

15.16±0.08 

3.65±0.10 

SOt, V2 = 1 

144,10 - 143,11 

360133.2 

898.7 

1.43 

X 

10 -“ 

4.42±0.37 

11.53±0.58 

13.44±1.23 

0.31±0.06 

so. 

345,29 - 344,30 

360290.4 

612.0 

4.80 

X 

10 -“ 

18.63±0.32 

11.14±0.11 

13.93±0.28 

1.26±0.06 

s‘«o 

IO 9 - 98 

360637.9 

90.6 

6.10 

X 

10 -“ 

11.30±0.52 

9.65±0.32 

13.98±0.71 

0.76±0.09 

S 02 * 

208,12 - 217,15 

360721.8 

349.8 

7.72 

X 

10-5 

15.07±0.55 

11.51 ±0.22 

12.59±0.56 

1 .12±0.12 

CH 30 H* 

11o,ii-10i,9£; 

360848.9 

166.0 

1.21 

X 

10 -“ 

8.12±0.44 

8.42±0.42 

15.26±0.87 

0.50±0.08 

NO 

/=7/2-5/2 

360948.3 

209.4 

4.52 

X 

10 -® 

15.30±0.55 

11.31±0.52 

29.18±1.19 

0.49±0.06 

CH 30 H* 

81 -72 E 

361852.3 

104.6 

7.70 

X 

10-5 

11.25±0.38 

6.58±0.27 

16.60±0.68 

0.64±0.07 

DCN* 

5-4 

362045.8 

52.1 

2.25 

X 

10-5 

2.26±0.33 

6.99±0.75 

9.17±1.22 

0.23±0.09 

3^802 * 

6^,3 “ ^2,4 

362158.2 

40.7 

3.29 

X 

10 -“ 

19.93±0.27 

11.17±0.08 

12 .88±0.21 

1.45±0.09 

”802 

152,14 - 14i,13 

362487.6 

120.8 

2.98 

X 

10 -“ 

6 .11 ±0.25 

10.16±0.26 

13.16±0.61 

0.44±0.04 

HNC,V2=1 

4-3 

362554.4 

709.3 

2.15 

X 

10-5 

4.20±0.29 

12.65±0.35 

11.39±1.06 

0.35±0.05 

HNC* 

4-3 

362630.3 

43.5 

2.30 

X 

10-3 

52.96±0.56 

5.72±0.06 

12.93±0.18 

3.85±0.42 

P-H 2 CO* 

5o,5 - 4o,4 

362736 

52.3 

1.37 

X 

10-5 

60.08±0.26 

6.91±0.03 

12.45±0.06 

4.53±0.25 

”80, 

232,22 “ 231,23 

362834.1 

259.2 

1.77 

X 

10 -“ 

6.75±0.31 

11.91±0.26 

11.59±0.65 

0.55±0.06 

SO, * 

214,18 “ 213,19 

363159.3 

252.1 

4.00 

X 

10 -“ 

50.66±0.33 

10.33±0.05 

15.12±0.12 

3.15±0.10 

”802 

234,20 - 233,21 

363286.4 

297.1 

4.01 

X 

10 -“ 

6.25±0.35 

9.08±0.35 

13.79±1.02 

0.43±0.07 

CH 3 OH* 

I 61 - I 60 

363440.4 

332.6 

2.40 

X 

10 -“ 

6.63±0.39 

9.74±0.46 

15.26±0.94 

0.41±0.06 

CH 3 OH* 

72 - 61 £ 

363739.8 

87.3 

1.70 

X 

10 -“ 

17.91±0.26 

8.36±0.10 

13.80±0.22 

1.22±0.07 

HC 3 N 

40-39 

363785.4 

358.0 

3.85 

X 

10-5 

3.67±0.22 

12.39±0.28 

9.71±0.70 

0.35±0.05 

SO, 

241,23 - 24o,24 

363890.9 

280.5 

1.78 

X 

10 -“ 

23.45±0.24 

10.42±0.07 

14.59±0.18 

1.51±0.07 

S 02 * 

232,22 “ 231,23 

363925.8 

259.9 

1.83 

X 

10 -“ 

27.90±0.34 

10.33±0.08 

15.27±0.22 

1.72±0.07 

P-H 2 C 0 * 

52,4 - 42,3 

363945.9 

99.5 

1.16 

X 

10-5 

25.65±0.28 

7.73±0.07 

12.59±0.15 

1.91 ±0.07 

P-H 2 C 0 

54,2 - 44,1 

364103.2 

240.7 

4.99 

X 

10 -“ 

11.15±0.28 

8.12±0.17 

13.31±0.39 

0.79±0.06 

0 -H 2 CO* 

53,3 - 43,2 

364275.1 

158.4 

8.88 

X 

10 -“ 

39.52±0.32 

7.79±0.07 

12.00 

3.09±0.14 

0 -H 2 CO* 

53,2 - 43,1 

364288.9 

158.4 

8.88 

X 

10 -“ 

39.70±0.38 

7.22±0.08 

12.00 

3.11±0.14 

Atmospheric 

HC3NV7=1 

7=40-39 

364403.7 

364676.3 

679.7 

3.86 

X 

10-5 

3.05±0.37 

15.90±1.10 

18.08±2.55 

0.16±0.06 

OCS* 

30-29 

364749 

271.4 

1.42 

X 

10 -“ 

8.45±0.31 

9.89±0.27 

14.07±0.56 

0.56±0.07 

H 3 O+ 

32,1 “ 22,0 

364797.4 

139.7 

2.78 

X 

10 -“ 

17.88±0.35 

8.58±0.14 

14.41±0.30 

1.17±0.09 

SO, 

259,17 “ 268,18 

364950.1 

497.1 

8.70 

X 

10-5 

10.11±0.26 

12.12±0.14 

11.08±0.34 

0.86±0.08 

CH 3 OH, v,=l 

53 - 52 £■ 

364986.8 

452.1 

3.76 

X 

10-5 

3.59±0.24 

6.02±0.36 

10.52±0.80 

0.32±0.06 

HNC,V2=1 

4-3 

365147.5 

709.6 

2.20 

X 

10-5 

3.79±0.60 

13.27±0.21 

7.70±1.42 

0.46±0.07 

HC 3 N V 7 =l 


365195.2 

680.2 

3.88 

X 

10-5 

3.00±0.28 

13.77±0.52 

10.33±1.05 

0.27±0.09 

P-H 2 CO* 

52,3 “ 42,2 

365363.4 

99.7 

1.18 

X 

10-5 

44.27±0.74 

7.73±0.12 

14.23±0.28 

2.92±0.16 

H 2 CN 

LO 

p 

1 

0 

365443.6 

52.7 

4.15 

X 

10-5 

2.98±0.31 

10.05±0.49 

9.13±1.07 

0.31±0.06 

H2C”S 

111 ,11 - lOi.io 

365613.4 

118.5 

7.33 

X 

10 -“ 

4.32±0.37 

9.25±0.51 

11.83±1.25 

0.34±0.07 

CH 3 COCH 3 

27ii,i6,0-26i2,15,1 

365696.2 

281.7 

1.20 

X 

10-5 

2.97±0.41 

9.46±0.91 

13.63±2.40 

0.20±0.06 

”S02 

241,23 - 24o,24 

365794.6 

279.9 

1.75 

X 

10 -“ 

4.56±0.45 

13.44±0.26 

7.54±0.68 

0.57±0.07 

U-line 

SO 2 , V 2 =l 

94,6 - 93,7 

365811.5 

365904.5 

842.5 

3.39 

X 

10 -“ 

6.31±0.84 

3.42±0.53 

10.91±1.14 

11.84±0.91 

19.03±3.20 

11.60±2.68 

0.31±0.07 

0.28±0.06 

SO 2 , V 2 =l‘’“/' 

74,4 - 73,5 

366125.8 

826.6 

3.05 

X 

10 -“ 

3.04±0.36 

12.82±0.71 

10.00 

0.29±0.07 

S02,V2=1'’“/2 

174,14 - 173,15 

366145.1 

943.8 

3.95 

X 

10 -“ 

4.28±0.32 

14.10±0.50 

10.00 

0.40±0.07 

SO,, V2=l‘’“” 

64,2 - 63,3 

366159.5 

820.0 

2.76 

X 

10 -“ 

3.68±0.34 

13.16±0.67 

10.00 

0.35±0.07 

SO, * 

15a,i4 - 14 i,i3 

366214.5 

119.3 

3.04 

X 

10 -“ 

91.77±1.56 

9.43±0.13 

15.80±0.32 

5.46±0.15 

HfCO 

5i,4 “ 4i,3 

366270.2 

64.6 

1.36 

X 

10-3 

6.96±0.77 

10.14±1.06 

14.00 

0.47±0.08 

”S02 

63 ^-^ — 5'7,4 

366521.9 

41.6 

1.31 

X 

10-5 

7.52±0.33 

11.61±0.25 

11.79±0.64 

0.60±0.07 

H 2 CN 

52,4 “ 42,3 

366637 

100.4 

9.42 

X 

10-5 

2.19±0.32 

9.06±1.04 

10.00 

0.21±0.06 

H 2 CN 

52,4 “ 42,3 

366655.4 

100.4 

7.69 

X 

10-5 

2.96±0.33 

8.36±0.80 

10.00 

0.28±0.06 

Atmospheric 

”S02 

254,22 “ 253,23 

366790.0 

367369.3 

337.6 

4.14 

X 

10 -“ 

8.97±0.37 

11.97±0.26 

12.69±0.60 

0.66±0.09 










Continued on next page 


5 Probably blended with 3^S02 ISi.w - 14i,i3 at 358988 MHz 
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Table A.l. The detected species toward the central position of the SLS field. The transitions marked with * have been detected toward AFGL 
2591 in the SLS survey (Van der Wiel 2011). The transitions marked with * have a non-gaussian line profile. The FWHM line widths of these 
asymmetric lines are from a Gaussian fitting. Their intensity is obtained by integrating over the whole line and their Vlsr is the velocity at Tpeak- 
The parameters without error bars are fixed parameters in the fit. 


Molecule 

Transition 

Freq. 

(MHz) 

(K) 

^9 
(s ') 



J Linbdv 

(K km s-*) 

Llsr 
( km s-') 

FWHM 
(km s-*) 

Lpeak 

(K) 

o-HzS*’*™'*/’ * 

32.1 “ 3i,2 

369101.4 

154.5 

1.90 

X 

kP” 

145.64±0.92 

8.54±0.04 

15.67±0.12 

8.73±0.43 


43,1 - 42,2 

369126.9 

262.8 

2.20 

X 

lO-*' 

15.23±0.76 

9.69±0.24 

10.45±0.56 

1.37±0.43 

HfS 

HfS 

32.1 “ 3i,2 

369170.4 

154.4 

2.55 

X 

10-5 

7.31±1.32 

11.26±1.40 

16.12±3.71 

0.43±0.13 

32,1 - 3i.2 

369246.1 

154.3 

1.91 

X 

lo-"* 

8.92±0.63 

11.15±0.42 

11.54±1.13 

0.73±0.15 

'^CS* 

8-7 

369908.6 

79.9 

1.06 

X 

10-5 

11.87±0.45 

10.00±0.26 

14.02±0.59 

0.80±0.10 

so, * 

96,4 - 105,5 

370108.6 

129.7 

4.19 

X 

10-5 

28.15±0.43 

11.41 ±0.09 

12.53±0.23 

2.11±0.12 

U-line 

44 - 84 

370499.6 

370931.7 

33.4 

1.45 

X 

10-5 

4.27±0.58 

2.55±0.39 

10.94±1.41 

9.63±1.24 

22.64±4.0 

13.90±2.46 

0.18±0.06 

0.17±0.08 

Atmospheric 
SO, * 

^■^,3 “ ^2,4 

371036.0 

371172.5 

41.4 

3.55 

X 

lo-’* 

81.33±0.57 

9.89±0.05 

15.66±0.13 

4.88±0.13 

SO,, V2=l 

214,18 - 213,19 

371264.8 

1017.1 

4.24 

X 

10-* 

4.11 ±0.71 

11.59±1.25 

13.72±3.49 

0.28±0.08 

CH 3 OH* 

1 

0 

> 

■H 

371415.7 

372.4 

2.52 

X 

lo-'' 

5.04±0.41 

9.57±0.62 

14.45±1.14 

0.33±0.09 

”S02 

254,22 - 253,23 

371804.6 

343.0 

1.35 

X 

10-*^ 

1.89±0.34 

10.69±0.85 

8.79±1.51 

0.20±0.08 

0 -H 2 CS* 

111 ,11 - 10i,io 

371847.4 

120.3 

7.72 

X 

lo-'' 

10.58±0.28 

9.06±0.18 

13.18±0.36 

0.75±0.09 

s‘’o 

99 - Sg 

372113.6 

102.5 

6.33 

X 

10 -* 

3.06±0.48 

9.84±0.95 

11.42±1.78 

0.25±0.09 

HNCO 

17i,i7 - 16i,i6 

372221.0 

204.1 

7.46 

X 

lo-’' 

8.13±0.64 

11.96±0.35 

9.23±0.83 

0.83±0.08 

^^SOj 

137,7 - 146,8 

372279.8 

199.7 

6.00 

X 

10-5 

1.73±0.37 

11.74±0.50 

6.30±2.23 

0.26±0.05 

N 2 H+ * 

4-3 

372672.5 

44.7 

3.33 

X 

10-5 

54.14±0.33 

7.44±0.04 

12.82±0.09 

3.97±0.20 

U-line'’ 

SO'^'O 

22o,22 “ 211,21 

372869.7 

373110.1 

221.2 

6.69 

X 

10 -“ 

7.92±0.41 

4.18±0.70 

9.33±0.26 

13.21±0.90 

10.05±0.64 

10.93±2.26 

0.74±0.11 

0.36±0.09 

SO* 

44 - 34 

373344.2 

33.8 

1.51 

X 

10-5 

43.76±1.28 

9.78±0.20 

14.32±0.51 

2.87±0.27 

HNCO 

17o,17 - 16o,16 

373600.7 

161.4 

7.54 

X 

10-“* 

21.05±0.55 

10.73±0.15 

11.84±0.35 

1.67±0.13 


HC 3 N 41-40 is probable, but the line intensity is not consistent with its other three detected transitions. 
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Fig. A.l. The identified lines toward the (from bottom to top) Eastern tail. Northern clump. South-west clump regions, and the center 
of W49A in the frequency range between 330 and 335 GHz and 335 and 340 GHz. 
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Fig. A.2. The identified lines toward the (from bottom to top) Eastern tail. Northern clump. South-west clump regions, and the center 
of W49A in the frequency range between 340 and 345 GHz and 345 and 350 GHz. 
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Fig. A.3. The identified lines toward the (from bottom to top) Eastern tail. Northern clump. South-west clump regions, and the center 
of W49A in the frequency range between 350 and 355 GHz and 355 and 360 GHz. 



Fig.A.4. The identified lines toward the (from bottom to top) Eastern tail. Northern clump. South-west clump regions, and the center 
of W49A in the frequency range between 360 and 364 GHz and 364 and 368 GHz. 
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Fig.A.5. The identified lines toward the (from bottom to top) Eastern tail. Northern clump. South-west clump regions, and the center 
of of W49A in the frequency range between 369 GHz and 374 GHz. 










